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ABSTRACT 


The electric field induced absorption lines of H 2 and D 2 have been 
studied. The frequencies of four Q-branch lines of H 2 and five Q-branch 
lines of D 2 have been measured as a function of density, and their 
shifts were observed to be in the linear region. The individual slopes 
and extrapolated zero density frequency of each line is determined. 

The frequencies of the lines were determined to be independent of the 
electric field intensity within experimental error. The polarizability 
of H 2 has been measured using the integrated intensity of the Qj (0) and 
Si CD , H 2 absorption line. 

A highly automated technique for determining the response function 
of the spectrometer using digitally recorded data is presented. The 
same techniques are used for manipulating data which is to be deconvo- 
luted and converted to absorption coefficient curves. For the Q x (0) 
and Q L Cl) lines of H 2 the halfwidths were measured as a function of 
electric field intensity at constant pressure. The half widths were 
observed not to be a function of the field intensity within experimental 
error. The half widths of these same lines were measured at several 
densities and compared to previously measured widths. 

Technical and operational details of several pieces of equipment 
built for this experiment, and equipment built for the five-meter Littrow 
spectrometer used for the experiment, are described. These include the 
Stark cell, a high voltage power supply, two carbon rod furnaces and 
their power supplies, one of which has been automated. Also included 

iii 
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are a high pressure gas filling system, the Stark cell fore optics and 
flushing cover, For the spectrometer, a detector pumping system, light 
chopper, remotely operated beam switching mirror and detector switching 
mirror have been. built. The modifications of the spectrometer optics 
to accept the Stark cell are discussed. 

An attempt to study the electric field induced spectra of N 2 was 
made. The unresolved Q- branch was observed but was too weak for any 
quanitative measurements. An unsuccessful attempt to form HD through 
a catalitic reaction of H 2 and D 2 is discussed. 
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CHAPTER I 


INTRODUCTION 

Diatomic molecules such as H 2 , D 2 and N 2 provide good subjects 
for study since their simplicity provides an excellent proving 
ground for molecular models and theories. Unfortunately such 
homonuclear diatomic molecules are normally infrared inactive due 
to their symmetry which places the center of positive and negative 
charge at the same position. By means of an externally applied 
electric field it is possible to distort the symmetry of these 
molecules slightly, not enough to perturb the energy levels of the 
molecule to any observable degree, but enough to induce a slight 
separation of positive and negative charge and therefore a small 
dipole moment through which transitions can occur. This technique 
works also with molecules with extremely small electric dipoles 
such as HD for which the electric dipole spectrum is practically 
unobservable. 

In electric field induced work the sample gas is introduced 
into a high pressure cell containing a pair of condenser plates with 
a high potential between them and with highly reflective surfaces. 
Light is transmitted through the gas between the plates by multiple 
reflections from the condenser surfaces. The sample gas is used at 
high pressures both to increase the number of absorbers in the path 
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and also to increase the electric field which can be maintained between 
the plates without electrical breakdown. 

This work contains the theory of electric field induced spectra 
as applied to diatomic molecules. In Section 2.1 the energy levels for 
the diatomic molecule is an electric field are rederived. The selection 
rules for electric field induced spectra are rederived in Section 2.2. 
They are AV = 0,1,2 — and AJ = 0,±2 for the vibrational and rotational 
quantum numbers respectively. The rotational selection rules are 
different from those found ordinarily in electric dipole transitions 
for diatomic molecules, and are identical to Raman selection rules. 

This difference arises from the fact that an external electric field 
applied in a certain direction induces a double angular dependence in 
the dipole moment integrals calculated for the electric dipole. This 
is the same angular dependence found in the dipole moment integrals 
calculated for the Raman effect and of course leads to the same 
selection rules. The selection rules are derived as one of the steps 
used in the determination of the strength factors which are also calcu- 
lated in Section 2.2. Section 2.4 calculates the Boltzmann weighting 
factors which are necessary to the determination of the integrated 
intensity. The expression for the integrated intensity contains the 
nonzero polarizability matrix elements through the strength factors so 
the integrated absorption is used to calculate the polarizability of 
the H 2 molecule as discussed in Section 4.6. The intensity of a line 
is also dependent on the orientation of the polarization vector of 
the incident radiation with respect to the direction of the electric 
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through the calculated strength factors, so that a polarizer must be 
used to eliminate one of the polarization directions when intensity 
measurements are being made. 

Since the gas is used at high pressures a shift in the line 
postions is introduced. For this reason the absorption line positions 
should be measured at several pressures whenever possible in order to 
extrapolate their positions back to zero pressure. This was done for 
the Q-branch lines of both D 2 and H 2 as discussed in Section 4,3 and 
4.4. A linear extrapolation was sufficiently accurate for the pressure 
range used in this experiment. 

The absorption lines are broadened by the high pressures used. 

For light molecules such as H 2 and D 2 this broadening is not as great 
as for heavier molecules such as N 2 . Section 4.7 and subsequent sections 
discuss the experimental procedure for the measurement of line widths 
including the determination of the response function of the spectrometer 
and the deconvolution process for removing the influence of the spec- 
trometer on the spectral line profile. 

The dipole moment induced in the molecule is proportional to the 
applied electric field intensity. Since the intensity of a line is 
proportional to the square of the dipole moment, the line intensity is 
proportional to the electric field intensity squared. High electric 
fields are therefore important and for intensity work, the field must 
be measured as accurately as possible. Chapter III contains a thorough 
description of the Stark cell and all associated equipment used in the 
experiment . 



CHAPTER II 


THEORY OF ELECTRIC FIELD INDUCED SPECTRA 

2. 1 Energy Levels of Diatomic Molecules in an Electric Field 

The energy levels of diatomic molecules are to a first approximation, 
the sum of the vibrational and rotational energies of the molecule. 
Transitions which take place between quantized energy levels are 
characterized by a change in the vibrational quantum number v and the 
rotational quantum number J. The spectral absorption lines resulting 
from these transitions consist of a series of rotational lines for each 
vibrational transition. It is found that the rotational energy depends 
on the vibrational state so according to Dunham (1), the molecular 
energy T can be expressed as 

T(v,J) = G(v) + F y (J) . 

where G is the vibrational energy and F the rotational energy. In order 
to obtain the energy levels, it is first necessary to derive the wave 
equation for the vibrating-rotating molecule. If we use the Born- 
Oppenheimer approximation we can separate the energies of the electrons 
from the energies of the nuclei. In the following only the nuclei are 
considered. 

If we place a diatomic molecule in a coordinate system fixed at 
the center of mass of the molecule whose component nuclei are m 1 and m 2 , 
as shown in the fbllowing figure, then the kinetic, energy is given quite 
simply by equation (1) . 
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2T = m r 2 + m r 2 . (1) 

112 2 

The center of mass condition and the definition of the vector r can be 
used to derive a one body problem from the given two body problem. They 
are 


m i*i + 


m 2 r 2 = 0 and 



From the equations above we can show that 


-> 

r 

1 



r and 


r 

2 




+in 


r. 


If we now substitute these equations into equation (1), the following 
one body problem results. 
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2T = 


m m 
—L-l 


m, +m„ 
l 2 


r z , 


If we define the reduced mass p as 


u = 


m i ni 2 

m 1 +m 2 


then 


T = 



The Hamiltonian for the system is then 


P 2 

H - * V « • 

where V(r) is the potential in which the nuclei move. Therefore 
Schroedinger' s time independent equation is 

- Iji V 2 U(r,0,<f>) + V (r)u(r,9,(j)) = EU(r,0,<j>) . (2) 

We can assume as usual that the wave equation U(r,0,0) is approximately 
the product of three wave functions, each depending on one variable only, 
so that 

U(r,e,<tO = R(r)0 ( 0 ) $ (<j>) . (3) 

Combining equations (2) and (3) and using the form of V 2 in spherical 


coordinates then 
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•h 2 

L- 1_ 

fr 2 3 1 . 1 

3 

r ■> 

ci nfi ^ 

i 

»*1 

2 y 

r2 3r 

3rJ r 2 sin 2 0 

36 

0 30 

r 2 sin 2 e 

3 4>2 


+ V'R0$ = ER04> . 

y\7. C "1 T‘1 ^ A 

I£ we multiply by — ~ W e find the <f> dependent part of the equation 

separates and yields with a separation constant of -m 2 , the usual 
azimuthal wave equation 

*(♦) = e im4> . (4) 


We are left with 


2 M 


sin 2 0 3 

f 2 3 I 

R 3r 

r 9r J 


Sine a 
0 36 



— 

sine fy 

© - m 2 

1 

w * 



+ Vr 2 sin 2 0 = Er 2 sin 2 6 


If we divide by sin 2 9 the e dependent part of the equation separates 
and yields with a separation constant of -J(J+1), the associated 
Legendre equation. The solution to this equation is the associated 
Legendre polynomials, usually written as 


0 (e) = pj(cose) . 


The remaining part of the wave equation depends on R only and is 

•1) 8ir 2 y 


r2 dr 


Id 2 dRi 

— -r~\ + 


dr 


‘ -J(J+] 


+ - V(r) 


R = 0 . 


(5) 


( 6 ) 
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The details of the solution of this equation are presented by Pauling 
and Wilson (2). The most important parts of their treatment are outlined 
in the following. 

Since the potential is appreciably anharmonic, we can use the Morse 
potential given by 


V (r) 


D ; 



e -a(r-r e ) j 2 


where D is the disassociation energy and a is a constant characteristic 
of the molecule in question. The constant is an equilibrium inter- 
nuclear distance corresponding to the minimum of the potential function. 
The Morse potential of a typical diatomic molecule is shown in Figure 1. 
We would like to cast equation (6) in the form of the radial equation 
of the hydrogen atom whose solution we already know. This equation is 


L" + | 2 ^ +1 ^ - ljl/+ ( y - 1 - ) j L a 0 , 


where A-£-l must be an integer value in order for solutions to exist. 
To this end we let 

R(r) = 7 S(r) , 


. e " a (r-r e ) and 


A - J(J*l)h 2 

8ir 2 yre 2 


( 7 ) 
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Equation (6) becomes 


d 2 S 8 tt 2 ii 
dr 2 + h^~ 


■ 

r 

2 

-A 

6 

. r . 

+ E - D(l-x) 2 


S = 0 


( 8 ) 


We can finish changing the variable in equation (8) to x by noting that 


d 2 S ? o d 2 S 2 ds , 
= i — y + a z x 1 — and 


dr 2 


dx 2 


dx 


1 - 


lnx. 


ar 


-2 


( 9 ) 


The latter equation can be expanded in a Taylor series in x-1. The 
first three terms are 


f -s 

r 

e 


1 + -=- (x-1) + 
ar v J 
e 


3 

a 2 r 2 


1 

ar 


(x-l) : 


Therefore 


= c 0 + c ^ X + C 2 x 2 , 


( 10 ) 


where 


c o - A 


, 3 3 

1 - + — n t 

ar a z r z 
e e 


(id 


C 1 = A 


ar 


a 2 r 2 


and 


( 12 ) 
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Now substituting equation (9) and (10) into equation (8) we obtain 


0 , fE-D-Cn 2D-C 1 

d 2 S 1 dS 8tt 2 u . ° 1 


dx2 + x dx a^h 


- D - C 2 S = 0 . 


Comparing the hydrogen atom equation to our equation above we let 

S(x) = e“ ay y^ F(y) and (15) 

y = y x • (1^) 

where a, 8 and y are arbitrary constants which can be manipulated to 
change the form of (14). We can assign the following constants to 
terms in equation (14) to make the choice of the values of ct, 8 and 
Y easier in later steps. 



8tt 2 u ,-or. r \ 
e = ( 2D ~ C l) 


When equations (15) through (19) are put into (14) the result is 


d 2 F dF 2B+ll „ . r, ( ^-b 2 ] _ ' fe/yz.?. qBV 

aF * 37 bd' 2 “ J l r J 
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+ (a 2 - d 2 /y 2 ) J = 0 . (20) 

Comparing this equation to the hydrogen atom wave equation, we can make 
the following assignments to a, f} and y to change the arguments in 
equation (20) so that they have the same form. 

6 = b; y = d/a; a = 1/2 . (21) 


Equation (20) becomes 

t22) 

where 

v = - \ (f+1) and (23) 

f = 2b . (24) 


In order for a solution to (22) to exist it is necessary to have v an 
integer as is the case of the hydrogen atom. To solve for the energy 
we must make use of equations (7), (11), (12), (13), (17), (18), (19), 
(23) and (24). We can solve (17) for E and substitute in (18), (19), 
(23) and (24) . After simplification the result is 


E = D + C. 




c ' 

-1/2 




v + 2 


a 2 h 2 

8ir z u 


1 

+ 2 



13 


The terms involving C^/D and C 2 /D can be expanded keeping only squared 
terms or less. The equation for E then becomes 


E - C 0 ♦ C 1 * c 2 - D 


'1 

4 


) 


+ c c + c 2 

12 2 


ah /2D 
2tt / pi 


Id ' c i * V * 4i? 


3C r 

C 1 C 2 + -f 


2"1 


1 

V + j 


a 2 h 2 

8ir 2 y 


v + 2 


After putting in (11), (12) and (13) and dividing by he then 

A 2 


__ a /2D __ 

he ~ he Da 2 r 2 hc + 2irC / y [ V 2 


a /2D 
2trc / y 


1 1 ] 

3A 

L + 1 

a /2D SA 2 

ar a 2 r 2 

e e J 

2D 

i V 2, 

2t re / ]j 8D z a 2 r 2 
6 


a 2 h f + ll 2 

I v + ij • 


If we put in (3) and let 


a /2D , 

to = ■= — / — and 
e 2ttc r y 


hw c 
e 

X e ~ 4D ’ 


then 


E_ 

he 


= GO 


V + 


-wx 

e e 


1 

+ 2 


)2 


8t r 


iyr e^ 


J(J+1) 


H ] CM 
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- (0 


e 


f 

1 

'i 

1 

3h 2 L , 1 

ar 

e 

a 2 r 2 

e J 

16Ti z yr e 2 D ^ ' 2 y 


J(J+1) 


. h 3 

64Tr 4 u2r^ B Da y -c 


J 2 (J+1) 2 ♦ 


If we let 


r - 11 . 

e Stt^ I c ’ 
e 

= 3h 2 1 1 

a e 16 tt 2 I D ar “ ; 

e e e 


^e 128TT b u> 2 I 3 c 3 ’ 
e e 

5w h 4 

e e = 512 tt 4 I 3 Da z and 
e 


I = 
e 


ur e 


2 . 


then 


E_ 

he 



ll 

. l] 

<D 

3 

ii 

[ v + 2j - w e x e 

v + 2 


+ 6 e J(J + D - 


i 'i i) 

v + 4 JCJ+l) - D J 2 (J+1) 2 - V + f J 2 (J+1) 2 
e e l 


Collecting terms we find that 


E_ 

he 



y- N 

1 


i 1 

2 1 

f \ 

, i 

dll 

e 

i v H 

- o) x 
1 e A e 

r ?J 

+ oj y 
e' e | 

v + 2 
^ / 


+ — 



2 


+ J(J+1) 



a 


e 




+ 



‘ J 2 W*» 2 [ D e * 6 e [ v *l) 

+ J 3 (J+1) 3 Jh 6 + — ] , 


+ 5 


l ] 2 

+ 2j 


+ 


where some higher order terms not kept in previous expansions have been 
added. If we define 


B = B - a 
v e e 


D = D + 8 Q 
v e e 


v + Ij + Y e 


M 


+ (5 


e l 


v + { 


+ — 


+ — 


H = H + — - and 
v e 


V V • 

^ = OJ 

r i 1 

v + T 

, \ i] 

2 

+ ui y 

f. ll 

V + 7T 

V->V 

- e 

2 

e e ( 2j 

e 7 e 

2 


where the subscript v denotes the vibrational dependence of the rotational 
constant, then 


|- = v > + B J(J+1) - D J 2 (J+1) 2 + H J 3 (J+1) 3 . (25) 

he v^v v v v v ' 

This expression yields the term values for the diatomic energy levels. 

If the allowed rotational and vibrational selection rules are substituted 
into the equation above, we obtain the allowed absorption transitions 
for the rotational branches within a vibrational transition. The 
selection rules for electric field induced spectra are derived in 
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Section 2.2 and are 


Av = 0, 1, 2 -- and AJ = 0,±2 , (26) 

where Av and AJ are changes in the vibrational and rotational quantum 
number respectively. The AJ = 0 transitions correspond to the Q branch; 
the AJ = +2 to the S branch and AJ = -2 to the 0 branch. The general 
formula for the frequencies in units of cm -1 of the three branches are 
therefore, combining equations (25) and (26) 

V J ) = Vv' + <V’V JCJ+1) " CV"V J2CJ+1)2 

+ (H v .-H v )J 3 (J+1) 3 , 

0 V (J) = V y ^ y , + B v .(J-2)(J-l) - B y J (J+l) - D y .(J-2) 2 (J-l) 2 
♦ D y J 2 (J + l) 2 + H y ,(J-2) 3 (J-l) 3 - H y J 3 ( J+ 1 ) 3 , 

S y (J) = + B y ,(J+2)(J+3) - B y J (J+l) - D y ^ (J+2) 2 (J+3) 2 

+ D y J 2 (J+l) 2 + H y ^ (J+2) 3 (J+3) 3 - H y J 3 (J+l) 3 . 

Here v is the initial (lower) vibrational state and v' is the final 
(upper) vibrational state. The J values are initial rotational states. 
Table I shows the first several transitions for specific J values, for 
the pure rotational spectral lines (Av=0) and the fundamental spectrum 
(Av=l) . The transitions in each branch are represented by a symbol 
which characterizes that specific transition. For example S 2 (3) 
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TABLE I 

ELECTRIC FIELD INDUCED TRANSITIONS 


Pure Rotational Transitions 


So(0) = 6Bo - 36Do + 216Ho 
So(l) = lOBo - 140Do + l,720Ho 

So (2) * 14Bo - 364Do + 7,784 Ho 

So (3) = 18Bo - 756 Do + 25,272Ho 
So (4) = 22Bo - l,364Do + 66,088Ho 

So(5) = 26Bo - 2,232Do + 148,616Ho 

So (6) = 30Bo - 5, 184Do + 299, 160Ho 


Fundamental Transitions 


QjCO) 

= 

V 0rH 



Qjd) 

= 

V 0->-l 

+ 

2B 1 

Q X C2) 

“ 

V 0-*l 

+ 

6B : 

Q l (3) 

= 

V 0+l 

+ 

12B 

QiW 

= 

V 0+l 

+ 

20B 

QjCS) 

= 

V 0+l 

+ 

SOB 

Q x C6) 

— 

V 0-*l 

+ 

42B 

SjCO) 


V 0+l 

+ 

6Bi 

s^D 

- 

V 0-»l 

+ 

12B 

S X C2) 

= 

V 0+l 

+ 

20B 

S l (3) 

=5 

v 0^1 

+ 

30B 

S l (4) 

= 

V 0t>-1 

+ 

42B 

SjCS) 

= 

V 0+l 

+ 

56B 

s^e) 

= 

V Q+1 

+ 

72B 


- 2 Bo - 4D X + 4Do + 8H Z - 8Ho 

- 6 Bo - 36D L + 36Do + 216H 1 - 216Ho 

l - 12Bo - 144Dj + 144Do + 1,728H 1 - l,728Ho 

1 - 20Bo - 400D 1 + 400Do + 8,000^ - 8,000Ho 

1 - 30Bo - 900D 1 + 900Do + 27,000H 1 - 27,000Ho 

1 - 42Bo - 1,764D 1 + l,764Do + 74,088^ - 74,088Ho 

- 36D 1 + 216H 1 

l - 2Bo - 144D L + 4Do + 1,728H 1 - 8Ho 
x - 6Bo - 400D 1 + 36Do + 8,OOOH 1 - 216Ho 
i - 12Bo - 900DJ + 144DO + 27,000H 1 - l,728Ho 

i - 20Bo - 1,764D 1 + 400Do + 74,088H 1 - 8,000Ho 

1 - 30Bo - 3,136D i + 900Do + 175,616H 1 - 27,000Ho 

- 42Bo - 5,184D 1 + l,764Do + 373,248H L - 74,088Ho 
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TABLE! (continued) 


°1 

(2) = v 

0+1 

- 

6Bo + 

36Do - 

216Ho 


Oi 

(3) = v 

0 + 1 

+ 

2Bi - 

12Bo - 

4Dj + 

144Do + 8H 3 - 1.728HO 

Oi 

(4) = v 

O’*' 1 

+ 

6Bj - 

20Bo - 

36D 1 + 

400Do + 216H 3 - 8,000Ho 

°i 

(5) = v 

0+1 

+ 

12Bj 

- 30Bo 

- 144D 3 

+ 900Do + 1,728H 1 - 27,000Ho 

Oi 

(6) = v 

0+1 

+ 

20B 1 

- 42Bo 

- 400Dj 

+ l,764Do + 8,000H 1 - 74,088Ho 
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represents a transition to the second vibrational state, the subscript 
2 giving the final vibrational state. Since in the transitions 
considered here the molecules are all initially in the ground state, 
the initial v state is understood to be 0 and is not specified in the 
transition symbol. The 3 in parenthesis represents the initial J state 
and the S implies a selection rule of AJ = +2. Therefore a rotational 
transition from J - 3 in the initial vibrational state 0 to J = 5 in 
the final vibrational state 2 is specified. 

2.2 Strength Factors of a Diatomic Molecule in an Electric Field 
The transition probability between two rotational states is 
proportional to the square of the dipole moment integral linking the 
two states. Since the intensity of an absorption line is directly 
proportional to the transition probability the squared dipole moment 
integral is referred to as the strength factor and is defined by 
Herzberg (3) as: 

S i (J,m;J"m") = | D. (J f m;J' f m') I 2 » (27) 

where D i is the dipole moment integral given by 

D i (J,m;J'‘,m") = j ^ RR \i i ip RR dx; i=x,y,z . (28) 

if; is the wave function of the rigid rotor and y . are the components 
RR i 

of the dipole moment in the x, y or z direction. Since the diatomic 
molecule is not strictly a rigid rotor equation (28) is an approximation 
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to the true equation. Before calculating the strength factors, 
it is necessary to calculate the dipole moment components in equation 
(28) for a diatomic molecule in an external electric field. 

According to Herzberg (3), when a molecule is subjected to an 
external electric field, a dipole moment is induced which is proportional 
to the field intensity. In laboratory coordinates where measurements 
are made the dipole moment is 


where i,j = ?c,y, z and a^. is the polarizability tensor. x,y,z are 
laboratory or space fixed coordinates assigned as shown in the following 
figure. 
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where E and E are zero in our assignment of coordinates. We find 
x y 


then that 




21 


y 

u 

y 


X 

y 

z 


- a E 
xz z 


* a E 
yz z 


= a E 
zz z 


and 


(29) 


The polarizability of a molecule is expressed in molecule fixed 
coordinates so we must transform the space components of the polariz- 
ability matrix to molecular coordinates. Let x' ,y' ,z' be assigned to 
the diatomic molecule as shown below, 


■y' 



with z" being the symmetry axis. Obviously the x" and y' axes are 
interchangeable. The transformation of the polarizability tensor is 
given by 


a 


ij 



A. 

is jt 


a s't' 


(30) 


The A's are the direction cosines describing the two rotations which 
are necessary for the transformation. According to Goldstein (4), A 
is given by 
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x 

A = y 
z 


x' y' z' 

cos<p -cos0 sin<)> sin6 sin<f> 
sirujp cos6 cos^ -sin0 cosif> 

0 sin0 cos0 

* ' 


(31) 


If we expand equation (30) and let the off diagonal terms go to zero, 
since the molecular axis is a principle axis system, then we find for 
the space fixed component a 

X. z • 


a = A ,A + A ^A + A „A . 

xz xx zx xx xy zy y y xz zz z z 


Using equation (31) , 


a s; (a * * - a „ ..) sin© cos0 sin$ 
xz z z yy 


Since x' and y' are interchangable by symmetry let us refer to these 
axes as IS to designate they are perpendicular to the molecular axis. 
Likewise we let z " be referred to as || Therefore 


a = (. a I, ' ~ a ) sin© cose sin<j> . 

X Z 1 1 ) J J- -L 

Similar expressions for the other two components of interest are 
obtained and our dipole moment components in equation (29) become 

u x = w' " sin0 cos9 sin( ^ E » 

p = - (ot ^ ^ - a,..,^) sin0 cos0 cos<t E and 

~y 1 1 -i- -1— 

P z = (a x ^a J _..sin 2 0 + cos 2 0) E . 
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If we introduce the definitions given by Brannon (5) for the anisotropy 
of polarization, y, and the average polarization, a, which are 

y = (a irjr - and 

a = j (2® x 'j_' * a n ',j > 

then we find 


y x = ysine cos6 sin<j> E, 
uy = -ysin0 cos6 cos<j> E and 
y z = [i Y (3 cos 2 0 - 1) + ct 


E . 


(32) 


These are the final expressions for the dipole moment components which 
are needed to calculate the strength factors as outlined in the following 
treatment . 

Manneback (6 and 7) has evaluated the strength factors for Raman 
spectra which also applies to electric field induced spectra. His 
treatment however ig more general than is needed for the specific -case 
at hand and some details are lacking. Therefore the strength factors 
were reevaluated for electric field induced spectra specifically and 
with a more straightforward treatment. We begin by evaluating the 

, y 

dipole moment integrals given by 


D i (J ,M; = 


*m “i *rr dr - 


where the rigid rotor functions are given by 


M 

1> RR = (cose.) *(40 . 
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Here 0 are the normalized associated Legendre polynomials and 4 is the 
usual azimuthal wave function. Since the space fixed coordinates are 
cylindrically symmetric, the x and y directions are indistinguishable. 

We can therefore let y^ = y ,| or y j_ where J_ refers to a direction perpen- 
dicular to the electric field and the )| direction is parallel to the 
electric field. Hence 


= Cm x 1 *V and 

• 

The dipole moment integrals for the parallel and perpendicular directions 
are, using equations (32), 

D (J ,M; J'jM/) = YEIxCeHxOfO 5 


lx(0) = 


0™ (9) sin6 cosQ 0'".. (0)dx (6) ; 
J J 






(<f>) [sin<fi + icosif>]$ m .((f')dT((fi) ; 


(33) 


D )5 (J NT) = [yElJ i;) (e) + (a-iy)Elf, 2) (6)]!,, (4) ; 


1^(0) = ©j"(9) cos 2 0 G m (0) dx (0) ; 

' J' 


ij 2 ) te) . 


e” (e) 0j,{0)dT(e) ; 


i„W) = 


4*(4>)4 m .(0)dT(<(!) 


(34) 
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Using 

*„,(♦) = 1/2 e1 ^ and dx(<i>) = d* , 

equations (33) and (34) can be evaluated. 


li ±3 (0) 

I, I (0) 


6 , . 
m,m 


These expressions give the selection rules for the azimuthal quantum 
number m, which are 


Am = 0,±1 . 


Our dipole moment integrals reduce to the following after substituting 
in the allowed m selection rules: 


D C i ) (J,mjJ f m + 1) = ^ yE Ix(6) ; (35) 


I±(6) = 9j (0) sinBcosfl (0)dT (6) ; (36) 


D (J,m; J',m) = yEl!, 13 (0) ♦ 


if. 1 '* (6) = 


a- jY 


Elf 2) (e);(37) 


Q™ (0)cos 2 0 0j.(0)dr(0); (38) 

if, 23 (6) = | ^*(e) ( e ) dT ( 6 ) . (39) 


Using the orthogonality condition for the associated Legendre polynomials 


(0)gP 1., (0)dx(6) = <5 T and dt(0) = sin6d0 
JO 0 J j J 


( 40 ) 
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equation (39) can be evaluated and yields 

1 n 23 C03 - « JfJ - • 

Equations (36) and (38) can be evaluated using recursion relations for 
the associated Legendre polynomials. These are given by Condon and 
Short ley (8) as 

rn nj HI 

COS9 Oj = A(J,m) 0 J+1 + B(J,m) 0 J _ J and 
m — 1714-1 m +1 

sine Gj = + C(J,m) Qj~\ ± D(J,m) ej^J , 
where the coefficients are given by 


Afj ml _ [(J-m+l) CJ+nfrl j1 1/2 . BfJ , „ r (J-m) (J»in) I 1 / 2 . 
A(J,n>) ( 2 J+1) (2J+3) J * B(J,mJ [(2J-1) (2J+1)J 


_ Rjim+l) (Jijn+2)! 1 / 2 . n , T _ R J+in ) (J+m-1)"] x / 2 

“ [_ (2J+1) (2J+3) J ' D(J,m) - [(2j_i) (2J+1)_ 


Using these recursion relations equation (36) becomes 

I x (9) = ”C(J',mTl)I(l) ± D(J', n*l)I(2) ; 

m* Tn 

1(1) = 0j cose 0j- +1 dx(8) ; 

j 

r 

K2) = ©J 1 * cose dt(e) . 


By reapplication of the appropriate recursion relation the two latter 
integrals break up into two integrals each 
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1 ( 1 ) = A(J'+l # m)I(3) + BGr+l,m)I(4); 

■(« • { ef e J' +2 dT < 9 ’ J 

1(4) - J 0j* Oj, dt(9) ; 

1(2) = A(J"-l,m)I(5) + B(J-l,m) ; 

1(5) - f Qj* Qj, dr(0) ; 

/ 

1(6) = f ©j* Bj ,_ 2 dx(0) . 

✓ 

The resulting integrals can be evaluated using the orthogonality 
relation (40) so that 

I x (0) = “ C(J',m7l) [A(J'+l,m)6 JjJ . + 2 + B(J'+l,iri) 

± D(J'’ > ni7l) [A(J"-l,m)6j y + B(J'-l,m) 6^ j,_ 2 ] . 

Using the same technique equation (38) can be evaluated with the result 
I ( ^(9) = A(J',m)[A(J"+l,m)6 J>J . +2 + B(J'+1 ,m) <5^ j,] 

+ B(J",m)[A(J'-l,m)6 JjJ , + B(J'-l,m)6 JjJ ._ 2 ] . 

Putting all the evaluated integrals back into the expressions (35) and 
(37), the dipole moment integrals become 

Y e |+C(J-\m+l) [A(J"+l,m)6 J ^ J . +2 

+B(J'+l,m)6 J y ] ± D(J',m7l)[A(J^l,m)6 J}J . (41) 

+ B(J'-l,m)<S J ^ J ._ 2 ]| ; 
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= y e|a(J% m) [A(J' + 1, m)<$ J><r+2 + 


+ B(J"+l,m) 5j J.] + B CJ ^ ,m) [A c j ' - 1 , nO s aJ ^ 

* + (a - |rj E Sj>J , . 


( 42 ) 


By inspection of (41) and (42) the selection rules for the rotational 
quantum number J are easily seen to be 


AJ = 0,±2 . 


If we substitute in the selection rules for all possible cases then 


D„ (J, m; J,ni) = yE[A(J, m)B(J+l, m) + B(J, m)A(J-l, m) ] + 


a - 3 Y 


E; 


D x (J, mjJ,m+l) = — Y E[C(J,m+l)B(J+l ,m) - D(J,m+l) A(J-l,m)] ; 


Dk (J, m; J+2,m) = yE B(J+2,m)B(J+l,m) ; 


Dj. (J,m;J+2,m+l) = ~j= yE 


D(J+2,m+l) B(J+1 ,m) 


(43) 


Du (J,m;J-2,m) = yE A(J-2,m) A(J-1 ,m) ; 

D x (J,m;J-2,m+l) = ^ yEC(J-2,m+l)A(J-l,m) . 

The strength factor. Si, is the square of these dipole moment integrals 
and is 


Si 


J.rnjJ'* (J) ,nT(jm) 


I = Di 


Di [J,m; J' (J) ,m" 



Since the m splitting is too small to be observed, we may average over 
them by summing over the m states and dividing by the degeneracy of the 
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particular J state. Therefore 


Si[J,J'(J)] = 


(2J+1) 


+J 

l 

m=-J 


Di[j, m ;J(0,m(m')] 


The summing of the expressions in equations (43) is quite tedious but 
straight forward. The final expressions are given in Table II. The 
sums involved in the evaluation of these expressions are given in 
Table III. The J dependent coefficients of the strength factors are 
listed in Table IV. 

2.3 Stark Splitting of m-Levels 

The energy of the split m-levels is, according to Herzberg (3), 

to = bi ■ -y p E , (44) 

o b 

where wo is the unsplit energy and y £ is the average component of the 
dipole in the direction of the electric field E. The electric field 
splitting is different from magnetic splitting in that equal values of 
m with different signs do not lead to different energies. Therefore a 
J=3 level would be split into 4 levels with 3 of them being doubly 
degenerate. From equation (44) the energy difference between the split 
m-levels is 


Aw = y E , 
z * 

where z is the direction of the electric field. From equation (32) of 
the previous section. 



TABLE II 


Sj.(J,J) = 
Sj,(J,J+2) 
Sj. (J >J + 2) 
Sj|(J »J-2) 
S ± (J,J-2) 


STRENGTH FACTORS 


a z E z i + 


4 J(J+l)y 2 E 2 
45 (2J-l)(2J+3) 


A_ Y 2 E 2 

15 Y (2J-1) (2J+3) 

- 2_ v 2p2 (J+l) (J+2) 
15 Y (2J+3) (2J+1) 

_ 1_ Y 2 E 2 (J+D (J+2) 

" 10 Y (2J+3) (2J+1) 


1 . Y 2 E 2 J(J-D 

15 Y E (2J-1) (2J+1) 


- i_ Y 2 £ 2 J(J-D - 

" 10 Y (2J-1) (2J+1) 





31 



h y / v p F (J-ro+i)(j+m 
s tj l Yb L (2J+1) (2J+ 


(J+ro+1) (J-m) (J+m) 1 F 

(2J+3) (2J-1) (2J+1) "3 


, r2J+ n r i_ * „2p2 ' 

+ 45 (2j_i)(2J+3) b 


H f l 1 / 2 

[ j (J±m) (J±m+1) (J-m+1) (J+m+1) (2J-1) 

«=-J ^ L J 


-|l/2 l 2 

(J±m+1) (Jim) (J-m) (J+m) (2J+3) \ 


= jj J(J+1) (2J+3) (2J-1) (2J+1) 3 


l (J-hh- 2) (J+m+2) (J-m+1) (J+m*l) 
m=-J 


= (2J+1) (J+l) (J+2) (2J+3) (2J+5) 


l (J+m+3) (J+m+2) (J-m+1) (J+m+1) = 
m=-J 


(2J+1) (J+l) (J+2) (2J+3) (2J+5) 


l (J-ta-1) (J+m-1) (J-m) (J+m) = yj J(J-l) (2J-1) (2J+1) (2J-3) 
m=-J 


l (J-iii) (J+m) (J±m-2) (J±m-1) = ~ J(J-l) (2J-3) (2J+1) (2J-1) 

m=-J 
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TABLE IV 

COEFFICIENTS OF THE STRENGTH FACTORS 


J 

J'=J 

Parallel 

J'=J+2 

Parallel 

J'=J-2 

Parallel 

0 

,000000 

.088889 

.000000 

1 

.035556 

.053333 

.000000 

2 

.025397 

.045714 

.017778 

3 

.023704 

.042328 

.022857 

4 

.023088 

.040404 

.025397 

5 

.022792 

.039161 

.026936 

6 

.022626 

.038291 

.027972 

7 

.022524 

.037647 

.028718 

8 

.022456 

.037152 

.029281 

9 

.022409 

.036759 

.029721 

10 

.022375 

.036439 

.030075 

J 

J'=J 

Perpendicular 

J'=J+2 

Perpendicular 

Perpendicular 

0 

.000000 

.066667 

.000000 

1 

.026667 

.040000 

.000000 

2 

.019048 

.034286 

.013333 

3 

.017778 

.031746 

.017143 

4 

.017316 

.030303 

.019048 

5 

.017094 

.029371 

.020202 

6 

.016970 

.028718 

.020980 

7 

.016893 

.028235 

.021539 

8 

.016842 

.028235 

.021961 

9 

.016807 

.027864 

.022291 

10 

.016781 

.027569 

.022556 
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u z = 


y y(3 cos 2 0 - 1) + a 


The average of u z is 


y 


z 


7T 

f 


y d6 

0^ 

z 

7T 

de 

0^ 




a 


E , 


so that 

E 2 . 

Using the polarizability of H 2 and a field of 250,000 volts/cm, this 
equation yields a splitting of about 0.0003 cm - *. This is of course 
too small to be detected, and does not even contribute to the width 
of the line (within experimental error) . It is therefore necessary 
to consider the m-levels as degenerate energy levels as was done in 
the previous section. 

2.4 Boltzmann Factors of a Diatomic Molecule 

The intensity of an absorption line depends on the number of 
molecules in the initial state, that is the number of molecules 
available for transitions to some upper state level. If we first 
consider vibrational population, the number of molecules in each 
vibrational state relative to the number in the ground state is 

\ _ -G 0 (v)hc/KT 
N 

o 


Ato 


~ 6 


Y + a 


( 45 ) 
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where N and N are the number of molecules in the vth level and the 
v o 

zero or ground state respectively. G q (v) are energy differences between 
the ground state and upper vibrational levels v. The total number of 
molecules is N and can be expressed simply as 


N=N q +N 1 +N 2 + , 


or using (45) 


N = N (1 + 
o v 


e -G 0 (l)hc/KT + e -G 0 (2)hc/KT + 


— ) • 


4 


If we let 


Q = 1 t e -Go(i:»>c/>CT * e -G 0 (2)hc/KT + ___ _ 


then 

N = N Q . (46) 

o 

Using (45) and (46) the total number of molecules in the vth state is 

N N e -G 0 (v)hc/KT _ 

V 

For the lighter diatomic molecules such as H 2 and D 2 the exponential 
quantity is of order lxlO" 9 so that for all practical purposes we can 
consider all molecules to be initially in the ground vibrational state. 


N = N . 
o 


so that 



The rotational population levels can be obtained by similar 
arguments if the degeneracy of the rotational levels is considered. 
Each rotational level has a degeneracy of 2J+1 which accounts for the 
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m states which are not resolved (see Section 2.3). Therefore 

N, . fi C2J + l)a- F t J ) hc/KT and 
J K 

R » l C2J,l)e- F t J > hc/lcr . 

J 


(47) 

(48) 


where 


F(J) = Bo J ( J + 1 ) - DoJ 2 (J+l) 2 + HoJ 3 (J+l) 3 . 

The ground state rotational constants are used since our molecule has 
been shown to always start in the ground vibrational state. An 
additional degeneracy exists in homonuclear diatomic molecules due to 
the identical nuclei. If I is the spin vector of the nuclei then T, 
the total nuclear spin, takes the possible values 


T = 21, 21-1, 21-2 --- 0 t 


For each possible total angular momentum there is a degeneracy of 2T+1. 
This degeneracy describes the number of wave functions necessary to 
describe the states. It is found that odd values of T and their 
associated degeneracy are associated with the odd J value wave functions 
and the even values of T with the even J value wave functions. This is 
due to the symmetry requirements of the total wave function. Therefore 
each odd and each even J state is subject to a further degeneracy 
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which we will call g(odd) and g(even). For example in D 2 the spin of 
each Deuterium nucleus is 1 . Therefore the total spin can be 

T * 2 or. 1 or 0. 

The odd T value 1 is associated with the odd J values and the degeneracy 
is 2T+1 = 3. The even T values 2 and 0 are associated with even J 
values and the degeneracy is 5+1 * 6. We therefore have a ratio of 
1 to 2 in the degeneracy of the odd to the even levels. Therefore for 

d 2 

g(odd) = 1 and 
g(even) = 2 . 

Likewise for H 2 where I = -j 

g(odd) = 3 and 
g(even) = 1 . 


Since for N 2> I = 1, the degeneracies are the same as for D 2 . 

Taking account of this degeneracy in equation (47) and (48) we 

have 

Nj - » and (. 

R = I g(2J*l)e- F t J > h ' :/l<T , 

J 


where 


g = g(odd) when J is odd, and 
g = g(even) when J is even. 
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Equation (49) can be written more conveniently as 

Nj = Nf j , (50) 

where f T the Boltzmann factor is 
J 

f _ g(2J+D -F (J)hc/KT 

Relative populations calcualted from equation (50) at T = SOO^K for 
H 2 , D 2 and N 2 are listed in Table V. From this table it can be seen 
that the population decreases very rapidly for increasing J values 
for H 2 and D 2 . For this reason relatively few lines in each rotational 
branch are seen. 

2.5 Integrated Absorption 

The exponential law of absorption for radiation passing through 
an absorbing medium of thickness l and with an absorption coefficient 
K is according to Herzberg (3) 



where I Q is the incident intensity and I is the intensity after 
absorption, 

In differential form 

dI ■ (S2) 

where N T is the number of molecules in the initial state and 

v,J v,J 

is the Einstein probability of absorption. From (51) and (52) we see 


that 
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TABLE V 

RELATIVE POPULATIONS 


J 

H 2 

°2 

N 2 

0 

.128695 

.181242 

.012623 

1 

.656396 

.204145 

.018578 

2 

.117728 

.384205 

.059616 

3 

.091875 

.114468 

.039419 

4 

.004275 

.094801 

.093943 

5 

.001014 

.014190 

.052206 

6 

.000016 

.006285 

.110099 

7 

.000001 

.000523 

.055608 

8 

— 

.000133 

.108272 

9 

— 

.000007 

.050996 

10 

_ _ _ 

.000001 

.093227 
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K = N t B V i J hv. 

V,J V,J 

This expression is valid for a single frequency spectral line centered 
at the frequency v- In actuality however we must replace K with an 
integrated absorption coefficient which is a function of frequency. 

We have therefore that 


a (v)dv = N t B V i J hv 
^ ' v,J V,J 


(53) 



The exponential law of absorption can be reapplied to the finite width 
lines with which we must deal and becomes 

-r r x , -a(v)& 

I fv) = I e v J or 
' o 

IfS- = , (54) 

0 

where we have assumed that I is independent of v. If the product 
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a(v)£ is small the right side of (54) can be expanded keeping the 
first two terms so that 

= i- a (v)£ . 

Solving for a(v) and integrating over the spectral line then 

| a(v)dv a £~ 1 j dv . (55) 

If the above approximation holds then it is said that the observed 
absorption is linear. Inspection of (55) shows that the observed 
spectral profile is directly proportional to the absorption coefficient 
curve if we disregard any instrumental distortion. If the product 
a(\>)£ is not small however then using (54) 

•In P~IcPj = “ a ^ £ ' 

Again solving for a(v) and integrating over v 
f ( ( ' 

a(v)dv = In dv , (56) 

The observed curve is no longer the absorption coefficient curve but is 
related through the logarithm of the fractional absorption. If we let 

r 

A = £ a(v)dv , 

where A is called the integrated absorption, then using (53) and (56) 
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A = 




J' 


hv 



Ihl] 

Io J 


dv . 


The integral on the right is measured from the observed curve; let us 
now expand the terms in the middle expression into physically observable 
parameters. is the number of molecules per unit volume in the 

initial state and was shown in Section 2.4 to be 





J 


where N is the total number of molecules per unit volume and fj is the 
Boltzmann factor also defined in Section 2.4. V* is the Einstein 

V j J 

probability of absorption and is given by 


8tt£ s v", J" 
v,J "he v,J 


Hor x . 


Here S 


v',J' 


v,J 


are the strength factors calculated in Section 2.2 

II or X 

which contain the squared dipole moment integrals and the electric field 


squared. If we let 
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A factor of three in the denominator is usually included in this 
expression in other experiments. This factor of 3: results from averaging 
over the polarization direction of the incident radiation. Since only 
one direction of polarization is used in this experiment for intensity 
measurements then 



orx 



NE 2 «j l 


v-,J- 

V,J 


V . 

fl orx 


The number of molecules per unit volume in the initial state is given 
by 


where n is the total number of molecules. The ideal gas law holds 
approximately for H 2 and other light molecules even for pressures as 
high as 450 p.s.i. so that 

N = £_ 

N KT * 


The final form of the integrated absorption is 


v',J' 8ir^ PE 2 A f 

v,J he KT J 

l| orx 


v',J 

V,J 


V 

ft orx 




(57) 


A is in units of cm" 1 . 

P is in units of dynes/cm 2 . 

E is in units of statvolts/cm. 


& is in units of cm. 
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^ is in units of cm 6 , 
v is in units of cm” 1 . 

T is in units of °K. 

2.6 Dicke-Doppler Width 

The Doppler contribution to the width of the spectral lines run 
in this experiment is according to a classical calculation 

Av = 2v | ~2KTln2~ l 1 / 2 
v c c |_ m J 

where Av c is the full Doppler width in cm 1 at half maximum intensity, 
v is the frequency of the spectral line also in cm m is the mass of 
the molecule and T is the absolute temperature. K and care constants 
and are Boltzmann's constant and the speed of light respectively. 
Classically therefore the Doppler width is independent of the gas 
pressure. A quantum mechanical derivation of the Doppler contribution 
to the width of a spectral line is quite different. For instance the 
derivations of Dicke ( 9 ) and Gersten and Foley (10) show that in the 
limit that the mean free path of the molecule is much less than the 
wave length of light, the Doppler width of a spectral line is 

Aw q = | t o -1 (£/X) 2 . (58) 

In this equation Aco q is the full width at half maximum intensity 
expressed as an angular frequency. t q is the mean time between collisions, 
T is the mean free path of the molecule between collisions and = A/2it 
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where X is the wave length of the spectral line. The condition that 
the mean free path is much less than the wave length is certainly met 
since for H 2 where X = 2.4xlO -4 cm(Q 1 (1) spectral line), the mean free 
path is about 3.8xl0" 6 cm at 20 atmospheres. If we express (58) in terms 
of wavenumber units and let t o _1 = f c , the collisional frequency, then 

. 4i:fcJl 2 v 2 

4v d 5T— ' 


where Av^is the full width at half maximum intensity and is called the 
Dicke-Doppler width. The collisional frequency is approximately 

fc % Nvo , 


where N is the number density, v is the average velocity of the molecules 
and a is the collisional cross section. The average velocity v is 


- 3KT 
v 

m 


1/2 


The mean free path can be expressed as the distance a molecule travels 
in time t divided by the number of collisions which occur in that time. 



_1 
No ' 


Therefore 


4irv 2 f3KT) 
Av D = 3cNcr m 


1/2 


In this calculation we see that the width of the line is inversely 
proportional to the pressure at constant temperature, since the number 
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density is proportional to the pressure. Hence the greater the pressure, 
the smaller the Doppler contribution to the line width will be. The 
Dicke-Doppler width calculated for the conditions and spectral lines 
used in this experiment are shown in Table VI. The Doppler contribution 
is seen to be small compared to the width of the lines encountered which 
were on the order of .050cm” 1 wide, and is therefore neglected in line 
width considerations. 

The calculated values are rough calculations using a cross 
sectional area derived from the hard sphere diameter which for H 2 is 
2.61 x 10” 8 cm. The hard sphere diameter of D 2 was assumed to be about 
the same as for H 2 . Since the Dicke-Doppler width is inversely propor- 
tional to the collisional cross section a larger cross section obtained 
by the inclusion of interacting fields would only decrease the width 
calculated here. The values calculated therefore comprise an upper 
limit to the expected Dicke-Doppler contribution. 
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TABLE VI 

DICKE-DOPPLER LINE WIDTHS AT 300°K 


Pressure 

(Amagats) 

Width of Q, (1) 1 
of H 2 (cm -1 0 

Width of Qj (l) 2 
of D 2 (cm -1 ) 

7.07 

.0052 

.0019 

10.12 

, 0035 

.0013 

13.17 

,0026 

.0010 

16.20 

.0021 

.0008 

19.22 

.0017 

.0006 

22.23 

.0015 

.0006 

25.23 

.0013 

.0005 

28.21 

.0012 

.0004 

31.19 

.0010 

.0004 


^Classical Doppler Width = .0363 cm 1 
Classical Doppler Width = .0185 cm -1 




CHAPTER III 


EXPERIMENTAL APPARATUS AND TECHNIQUES 

3.1 Stark Cell 

The Stark cell shown in Figures 2 and 3 contains the sample gas 
at high pressure between two parallel and highly reflective condenser 
plates. A large d.c. potential is maintained across the plates to 
provide a strong electric field in the sample gas. The condenser 
plates are ground steel held in a parallel plate arrangement by 
plexiglas holders. The plates are two inches wide and one-half inch 
thick with all edges rounded to minimize high local fields and the 

resulting corona and electrical breakdown. The plates are coated with 

/ 

a reflective nickel plating to transmit light as efficiently as possible 
by multiple reflections the length of the plates. The plates used in 
this experiment were 22-7/8" long and were held .067" apart by the 
plexiglas holder with no spacers between the plates. Many previous 
cells of this type used spacers between the plates. According to Boyd 
(11 and 12) these spacers limit the field which can be applied across 
the gas. The plates used in this experiment are the same that were 
used in previous work. The plates were badly pitted due to frequent 
breakdown and were re-ground and gold plated. 

The optical quality of the gold plating was far superior to the 
original nickel plating which was ground off but the strength of the 
plating was not equal to the task of withstanding the breakdown which 
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Figure 2. High Pressure Stark Cei! 
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either accidentally or purposely occurs. The gold plating was observed 
to be penetrated completely in several places by the breakdowns and left 
sharp edges which limited the field. The gold plating was removed and 
a new nickel plating was deposited which was of much better optical 
quality than the original nickel plating. As before "aging" the plates 
by breaking them down purposely appreciably increases the field which 
can be maintained between the plates. This process apparently burns off 
dust and high spots on the plates. The plates are held parallel by a 
plexiglas holder, and nylon screws thru the holders into the plates 
hold them in. place. The plate spacing was at one time increased to 
.075" to improve the transmission of the cell. This was decreased again 
due to the problems encountered: in the higher potential required to 
maintain the same electric field across the wider spacing. The plate 
spacing was measured by placing a small bit of wax between pieces of 
aluminum foil at several places on the plates which are then clamped 
together. The flattened pieces can then be measured with a micrometer 
to obtain the plate separation and its variation. The plate separation 
measured in this way was .067" ± .001". The entire plate assembly slides 
into the high pressure cell in a rectangular space created by cementing 
plexiglas fillers in the cell. These spacers also serve to reduce the 
unneeded volume of the cell. The cell is evacuated and filled through 
an inlet in the cell body. A valve is also fitted to the cell body 
to allow flushing. 

The cell body is made from thick wall steel pipe 3/4" thick. This 
is thicker than necessary to withstand the maximum 450 p.s.i. pressure 
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used but is used to provide adequate room for 0-ring grooves which seal 
the different lengths of the cell together. The cell is made in sections 
which can be interchanged to give different cell lengths. The thick wall 
pipe also provides adequate rigidity to the cell since it is bolted to 
the spectrometer at one end and must support its own weight. 

The end caps are 3/4" thick brass and contain the windows and high 
voltage connections to the plates. The windows were a continuing 
problem until the present solution was found. Originally very thick 
(5/8") windows were seated in a bed of silicone rubber in recesses in 
the end caps. It was found, however, that the windows bent under the 
high pressure and pulled away from the gasket material creating a leak. 
The end caps were changed to accept a much smaller and therefore thinner 
window sealed with an O-ring. These windows, however, broke when the 
cell was pressurized. It was discovered with the aid of a dial 
indicator that the center of the 5/8” thick, original end caps bowed 
.020" in the center under 400 p.s.i. pressure. This, of course, would 
bend the rigidly held windows until they broke. New end caps of 3/4" 
thick brass were made with window holders which "float" with the applied 
pressure. The window holder holds the window rigidly sealed against an 
0-ring and the holder itself is sealed against a second 0-ring. However 
as pressure is applied the window holder simply compresses this second 
0-ring more instead of bending with the end cap. No windows have been 
broken with the new holders shown in Figure 4. NaCl and CaF 2 windows 
were used and are 1/4" x 3/4" x 2". 
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The high voltage connector shown in Figure 5 passes through one 
end cap and is composed of a machines nylon body sealed with an O-ring. 
A brass feed-through also sealed with an O-ring carries the voltage 
into the cell. The connector inside the cell is rounded, as are all 
metal surfaces inside the cell, and plugs into a socket connected to 
one of the plates. The other plate is grounded to the cell body by a 
spring attached to the end cap. 

A thermistor was used to monitor the temperature inside the cell. 
This thermistor was soldered to the end of a feed-through shown in 
Figure 6. This feed-through was machined from a standard brass fitting 
and bakelite, and was sealed with epoxy. A bridge circuit with the 
thermistor as an active element was made and the current in this bridge 
network was calibrated against a mercury thermometer. The temperature 
in the cell was on the average 299. 8°K and did not vary more than 2°K. 
In calculations made in this work which involved the temperature, a 
rounded off value of 300°K was used. 

The entire cell is bolted to an access port on the spectrometer 
using a machined adaptor made for that purpose, and the other end 
containing the high voltage connection is attached to a flushing cover 
containing the source optics. 

3.2 High Voltage Power Supply 

The high voltage supply shown in Figure 7 is used to apply a large 
static potential across the condenser plates in the Stark cell. The 
field which can be applied is limited by the dielectric strength of the 
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compressed gas between the plates. The high voltage is supplied by a 
large d.c. transformer power supply from Plastic Capacitor. The voltage 
going to the cell is adjusted by means of a variac on the transformer 
primary, and is measured through a divider string with a ratio of 80,000 
to 1 ± 0.5% and read on a digital voltmeter. The majority of the 
resistors in the divider string are immersed in oil to prevent corona 
problems. About 300, 3.3 megohm, 1 watt carbon resistors were soldered 
together in series and strung between ceramic standoffs in a box filled 
with transformer oil as seen in Figure 8. The large number of resistors 
were used to prevent voltage drops across individual resistors which 
might lead to breakdown across the short physical length of the resistors. 
The small voltage drop also insures that the wattage rating of the 
resistors is never approached in the event of cell breakdown since this 
could change the calibration of the divider string. The voltage is 
read on a zero to one volt Simpson digital voltmeter. A 22 megohm, 
one watt resistor is placed in series with the high voltage line leading 
to the cell to limit the current to less than the five milliampere 
rating of the power supply in the case of electrical breakdown across 
the plates. This resistor and the connections leading to the cell and 
divider string are also immersed in transformer oil to prevent corona 
problems. Occasionally during breakdowns the wattage rating of the 
current limiting resistor was greatly exceeded causing the resistor to 
explode. In this way the resistor also acted as a fuse to protect the 
power supply. The connection to the cell is made with a brass plug 
rounded to prevent corona. The dry atmosphere inside the flushing cover 
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Figure 8. Divider, String 
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also helps with this problem. With the .067" plate spacing the 
operating- potential was about 30,000 volts, corresponding to a field 
of about 180,000 volts/cm. 

3,3 Gas Filling System 

The gas filling system shown in Figure 9, is designed to handle 
the high pressures used in this experiment and is also evacuable. The 
main manifold is brass, drilled and fitted with adaptors designed to 
accept large valves, good for pressure or vacuum. A Hastings Gauge 
thermocouple and mechanical pump are isolated by valves from the 
manifold when high pressure is being used. A high pressure bourdon 
type gauge calibrated in two-pound subdivisions can be isolated from 
the rest of the system if the manifold is being evacuated. Gas lines 

from the manifold. to the cell are 1/4" or 3/16" O.D. to withstand the 

high pressures. A liquid Nitrogen trap is included in the gas line to 
the cell. It is of course extremely important to get all the oxygen 
out of the manifold and cell before operating the reaction chamber 
or turning on the electric field, which might lead to an explosion. 

This is accomplished by evacuation and flushing. After the cell is 
filled, any condensable gases can be removed by the cold trap in the 

gas line to the cell if H 2 or D 2 is being used. It was found that N 2 

liquefies at 400 p.s.i. in a liquid Nitrogen trap which prevents its 
use in this case. 

The N 2 , H 2 and D 2 used in this experiment were from Matheson Gas 
Co., and the stated purity was 99.99% for N 2 and 99.9% for H 2 and D 2 . 




Figure 9. High Pressure Gas Filling System 







61 


A small secondary manifold was constructed to mix and D 2 and to 
contain a heated platinum wire to attempt to convert the mixture partially 
to HD. The mixing chamber consists simply of a common chamber into 
which H 2 and D 2 are admitted simultaneously from lecture bottles of 
equal volume and pressure. This is released at 400 p.s.i. into a 
reaction chamber shown in Figure 10. The reaction takes place via a 
heated platinum wire. The reaction chamber is stainless steel and is 
sealed with a teflon gasket. Teflon is used instead of an 0-ring since 
silicone grease used on an 0-ring will "poison" a platinum wire and 
render it almost inactive. Four feet of platinum wire .005" in diameter 
is wrapped around ceramic insulators and is connected to a variac through 
small glass feed-throughs. 15 volts at 1.1 amps are used to heat the 
wire. The gas is valved off in the reaction chamber and kept there for 
25 minutes while the platinum wire is heated, before being passed to 
the cell. This procedure is repeated until the cell' is pressurized. 

The platinum wire was obtained from the Poly Research Corporation. 

The equilibrium constants for the reaction : is according to Hill (13) 
given by the equation 

H 2 + D 2 t 2HD 

is according to Hill (13) about 3.5, so that if the reaction goes to 
equilibrium the final mixture of gases should contain almost 50 percent HD. 
At 445 p.s.i. total pressure this would correspond to a 215 p.s.i. 
partial pressure of HD which should be observable. Observation of the 
partial pressures of H 2 and D 2 in the gas mixture indicated approximately 
equal parts of each but far more than would be expected if the reaction 
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had gone to equilibrium. No HD absorption lines were observed. It is 
felt that the high pressures reduced the mean free path to the point 
that only relatively few molecules came in the vicinity of the platinum 
wire. 


3.4 Carbon Rod Furnaces 

Two carbon rod furnaces have been built for use as infrared 
radiation sources and are shown in Figures 11, 12 and 13. The furnaces 
consist mainly of a carbon rod held in large copper electrodes through 
which a very high current is conducted. The resistive heating of the 
rod causes it to reach a temperature of about 2500°K as measured with 
an optical pyrometer. This provides usable energy throughout the 
infrared and since the emmisivity of carbon is close to unity its 
radiation curve should closely resemble a black body curve. The rod 
is slotted in the center to increase the resistance of this portion of 
the rod and therefore increase the temperature. This slotted portion 
is the part of the rod which is used as a source. 

The carbon rods used in the source are spectroscopically pure 
carbon purchased from Ultra Carbon Inc. and are 5/16" in diameter and 
3-3/4" long in the first furnace built and 5/16 by 4-7/8" in the new 
automatic furnace. The rods are slotted in the center at a length of 
1-1/2" in the furnace used in this experiment and 2-1/2" in the other 
furnace which was made to bolt to the spectrometer. The cross section 
of the slot is V-shaped and is made with a 60° end mill cutter to a 
depth of .110". The V- cross section was found to have a longer life 
than a rectangular cross section slot also tried. Another advantage 
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Figure 11. Carbon Rod Furnace 
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Figure 12. Automated Carbon Rod Furnace Side View 
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of the V-slot is that a slight misalignment in the rotational placement 

/ 

of the rod in its holder should make no difference. 

i 

The rod is held in place by two, one-inch diameter copper electrodes 

i 

which carry power to the rod. Because of the high temperature of the 
rod the electrodes are water cooled by means of holes drilled in the 
electrodes. The entire body of the furnace is also water cooled. In 
the older furnace this is accomplished with copper tubing soft soldered 
to the outside of the case and to the top and bottom. In the newer 
furnace this is also done on the main body of the furnace but the front 
and back plates have been made hollow to allow water circulation for 
cooling. I 

The case is flushed with argon before turning on the furnace and 
kept at about 5 to 10 p.s.i. positive pressure (gauge) while the rod is 
operating. The gas line to the argon supply is left open to maintain 
the pressure against any leaks. Some form of drying agent is used in 
the gas supply line. Water inside the case drastically shortens rod 
life as was evidenced by small accidental cooling water leaks. The 
drying agent also protects the hydroscopic windows used in the furnaces. 
NaCl and KBr windows have been used and withstand the high temperature 
very well, as long as no large sudden temperature changes take place. 

3.5 Carbon Rod Furnace Power Supplies 

Because of the large power requirements 220 volt transformers are 

used in the carbon rod furnace power supplies shown in Figures 14 and 

I 

15. The secondary of the transformer will deliver up to 500 amps at 
12 volts. The primary voltage ia adjusted by a large variac rated at 
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Figure 14. Carbon Rod Furnace Power Supply 
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Figure 15. Automated Carbon Rod Furnace 
Power Supply , 
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25 amps. The current in the secondary is measured by a loop type 
current transformer and Simpson meter calibrated to read 500 amps full 
scale. The furnace is normally operated at 400 amps. 

A safety feature included in both power supplies automatically 
turns off the power to the rod in case the cooling water is turned off 
or drops significantly in pressure while the furnace is in operation. 

Also the rod cannot be heated without turning on the water. This is 
accomplished by passing the primary current through a large normally 
open relay before going to the transformer. The relay coil closing 
this relay switch can be energized only if a diaphragm type pressure 
switch is held closed by the cooling water pressure. The pressure 
switch is placed in the exhaust water to insure that a broken or leaking 
line anywhere in the series of cooling lines will open the relay and 
shut off power to the rod before any damage can occur to the furnace. 

To prevent excessive water from being spilled in case of a water line 
leak or breakage the pressure switch also' controls a solenoid in the 
water supply line which shuts off the water supply. A bypass around 
the pressure switch is necessary to initially set the system. 

The current to the rod is increased and decreased slowly by means 
of the variac to prevent thermal shock to either the carbon rod or the 
furnace windows. About 15 minutes is sufficient time to turn the 
furnace up to 400 amps. 

The newer power supply has been automated so as to free the operator 
from the tedious chore of turning the current up and down slowly. This 
was accomplished by attaching a properly geared down motor with a chain 
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drive to the variac shaft. Microswitches limit the variac travel to 
the proper current. So to turn the furnace on, it is necessary only 
to throw a switch which delivers power to the variac motor after the 
bypass button is depressed to initially set the pressure switch. To 
turn the furnace off the switch is reversed which reverses the variac 
motor. A timer mechanism is then automatically actuated which shuts 
off power to the relay coil after the variac motor has had sufficient 
time to turn the current down. Changing a rod in either furnace is a 
fairly rapid process. After the power is turned off the access flange 
is removed and four Allen screws are removed to replace the rod which 
is pressed into a copper collar at each end of the rod to provide good 
contact between the rod and the electrode. Periodic cleaning of the 
contact surfaces is necessary to maintain good contact. Poor contact 
results in a very eratic and lower rod current. The surfaces inside 
the furnace should also be wiped off periodically. These surfaces 
become coated with a carbon film and there is some evidence that this 
film conducts some current away from the insulated electrode to the 
grounded case. This was observed on one occasion on the newer furnace 
after it had not been cleaned for some time. The source was noticeably 
dimmer at the regular operating current. After cleaning off the accumu- 
lated carbon film, the original light intensity was recovered. This 
was confirmed by temperature measurements with an optical pyrometer. 
Often when a rod breaks arcing across the break will occur. This 
accelerates the deposit of the carbon film and can quickly cause a 
decrease in the intensity of the light if it is not removed. 
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3.6 Flushing Cover 

The flushing cover is a sheet metal box enclosing the source optics 
which is sealed and used to flush atmospheric gas from the source to the 
cell. The flushing cover is attached with machine screws to a table 
built to hold the source optics and is sealed with silicone rubber 
cement and Apeazion-Q. A sealable plexiglas access port allows working 
inside the flushing cover. An ordinary light bulb, switched on and off 
from outside the flushing cover, is kept inside near the hydroscopic 
window of the carbon rod furnace. This light serves to illuminate the 
interior of the cover for working inside and also keeps the temperature 
of the window above the ambient temperature of the surroundings and 
reduces condensation of atmospheric water on the window when the flushing 
cover is not purged. 

Nitrogen is used to flush the atmospheric gas from the flushing 
cover. The nitrogen is kept at a slight positive pressure to keep 
unwanted gases from re-entering the cover when 'it is t n6t purged. 

Nitrogen is used to flush the atmospheric gas from the flushing 
cover. The nitrogen is kept at a slight positive pressure to keep 
unwanted gases from re-entering the cover since it leaks slowly. The 
pressures used are up to 10 inches of water (less than 1/2 p.s.i.). 
Because this low pressure is difficult to regulate from a high pressure 
gas cylinder and since even a few pounds could exert several hundred 
pounds total force on the large surface area of the cover, a safety 
device was installed to limit the pressure to 10 inches Of water. 
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A solenoid operated valve is used in the gas line leading to the 
flushing cover. The solenoid is normally closed and is therefore "fail 
safe" in case of a power failure. A variable diaphragm type pressure 
is actuated by the pressure inside the cover and is set to open at 10" 
of water cutting off the power to the solenoid and shutting off the 
gas supply. When the pressure falls below 9" of water the gas line 
is re-opened. Since the solenoid opening and closing causes noise 
transients in the electronics, the pressure is kept below 9" of water 
so that the solenoid remains open, unless the pressure accidentally 
exceeds 10" of water. 

Several other devices are used in addition to flushing for removing 
water and carbon dioxide from the light path. Two open containers of 
sodium hydroxide and alumina are kept in the flushing cover. Sodium 
hydroxide absorbs water and carbon dioxide and alumina absorbs water. 

In addition to these chemicals, another device for removing condensable 
gases has been installed. A cold trap of 3/4" copper tubing with a 
valve at each end is attached to the flushing cover with one end opening 
through one valve to the interior of the flushing cover and the other 
opening through the other valve outside the cover. When the valve to 
the cover is opened and the trap immersed in liquid nitrogen, condensable 
gases are collected in the trap. If this valve is closed after a few 
minutes and the liquid nitrogen removed then the gases can be exhausted 
through the other valve. The gas removed is replaced with nitrogen. 

3.7 Spectrometer and Optics 

The spectrometer used in this experiment is a high resolution 
5-meter Littrow built at the University of Tennessee shown in schematic 
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in Figure 16, Details of the construction and operation of the 
spectrometer are given by Jennings (14) • The spectrometer is constructed 
in four separate optical sections, each evacuated. The sample section 
is designed to hold source optics and sample cells. The Stark cell 
used in this experiment is bolted to an access port cover outside the 
spectrometer by means of an adaptor plate made for that purpose. The 
cell is operated outside of the evacuated spectrometer since the high 
voltage, high pressure and normal maintenance make it impractical to 
operate the cell in a vacuum. The access port cover has been stiffened 
with steel members to prevent any flexing of the port cover and subse- 
quent movement of the cell bolted to it when the sample section is 
evacuated. O-rings seal the cell to the access port. 

The external cell necessitates an external source and source optics. 
The source used is a carbon rod furnace described in Section 3.5. The 
source is bolted and sealed with an O-ring to the flushing cover 
described in Section 3.6. The flushing cover contains the source optics 
and makes it possible to flush atmospheric gas from the path from the 
source to the cell which is also bolted and sealed with O-fings to the 
flushing cover. 

The source optics are shown in Figure 17. A cylindrical lens is 
placed in front of the source and is used to get more light through the 
Stark cell than would be possible without the lens. Theoretically 50% 
more light would be possible, in practice about 30% more was achieved 
in comparison tests. The cylindrical lens is Calcium Fluoride 1-1/2" 
X2-1/2" with a focal length of 3" and was obtained from John Unertl Inc. 





Figure 17 Source Optics and Stark Cell 
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The lens does not affect the vertical diverging light coming from the 
source, but reduces the divergence of the horizontal light. The 
horizontal diverging light striking the focusing mirror (Mj) is focused 
on the end of the cell, as is necessary in this type cell, but the 
vertical light is delayed in focusing until the center of the cell. 

In this way double the solid angle in the vertical direction can be 
accepted by the cell. 



The light from the focusing mirror is directed to the cell via a flat 
(M 2 ) which provides the proper distance to the cell for focusing. The 
light enters the spectrometer and passes through the optics (M 3 to M g ) 
designed to contain a sample cell. At a focus a mirror (M g ) has been 
installed on a raising and lowering apparatus. This mirrQr elevator 
was adapted from a cell elevator previously used in this laboratory. 

The mirror is attached to a plate riding in a dovetail track and is 
driven up and down by a rack and pinion driven by a servo. Microswitches 
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limit the travel of the mirror. The cycle time was decreased from 
about a minute to five seconds, and so to slow down and stop the mirror 
after the driving servo has been turned off by the limit microswitches, 
springs were installed which stop the mirror very effectively. The 
slave servo is driven by a master servo outside the spectrometer which 
is in turn driven by a geared down motor. Thus, only a switch need 
be activated to move the mirror up or down. The mirror enables either 
the primary or Stark cell beam to be used. With the mirror up, the 
light from the primary beam is reflected from the mirror and goes to 
the prism pre-disperser section. With the mirror down, the light 
from the Stark cell beam is allowed through and follows the same path 
through the spectrometer. In this way two experiments can remain set 
up in the same section. This setup also allows one section to contain 
an experiment and the other a calibration gas. It is for this purpose 
that the cycle time was decreased so that only a small amount of 
information is lost in switching from the experimental beam to the 
calibration beam. The primary beam contains a long path white cell 
being used in an experiment but was "borrowed" when calibration lines 
were needed. Carbon monoxide absorption lines were used for calibration 
with pressures of only 2 to 5 millimeters needed in the 12 meter path 
for the overtone band. 

The prism predisperser section contains the light chopper described 
in Section 3.8 and the optics necessary to double pass the prism. Both 
NaCl and LiF 2 prisms were used. From the prism predisperser the light 
is passed into the monochrometer and is focused onto the entrance slit. 
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The monochromator optics are designed to double pass a 10" x 16" 
grating. Both a 31.6 groove per millimeter and a 79 groove per 
millimeter grating blazed at 63° and 65° respectively, were used. 

After passing through the exit slit, the light strikes a mirror which 
has been mounted on a rotating table made from old spectrometer parts. 

A slave servo turns a worm gear which turns the table and is linked 
to an external hand turned servo allowing the table to rotate the mirror 
so that the light from the exit slit strikes either one of two detectors 
available in the detector section. Since the angle through which the 
detector switching mirror turns is very small, the driving servo is 
turned with a vernier type adjustment to allow fine angular adjustments. 

The detector section holds two detector cryostats and detectors. 
Liquid nitrogen cooled PbS and InSb detectors were used in this 
experiment. The detectors are from Santa Barbara Research and are 
contained in vacuum insulated dewars. The dewar pumping system is 
described in Section 3.9. 

3.8 Light Chopper 

A barrel type light chopper was designed and constructed and is 
shown in Figure 18 and 19 and used to chop the radiation reaching the 
detectors, and to produce a reference signal in synchronization with 
the radiation chopping frequency for the amplification electronics. 

The chopper barrel has 18 slots and is driven by a 3600 rpm synchronous 
motor to chop the radiation at 1080 Hertz. The slots are 2-1/2" tall 
to accommodate the tallest image the spectrometer is designed to accept 
in this section of the optics. The barrel diameter, the f-number of the 
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18-Tooth 

Gear 


Figure 19. Light Chopper Top View 
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beam focused at the center of the barrel, and the chopping frequency 
were all considered as variables in arrving at a barrel size and number 

i 4 

of slots. The barrel diameter is 2.865" and the slots are 1/4" wide 
with 1/4" between the slots to give an approximate sine wave output. 

The barrel and barrel shaft were made in one piece to assure good 
alignment of their axes and the barrel was made as thin as practical 
(.030") to reduce the inertial loading of the motor. A steel 18-tooth 
gear on the shaft with the chopper barrel turns between two iron pole 
pieces carrying the field from a permanent bar magnet. As the gear 
teeth pass between the poles a small emf is picked up and transformed 
by a coil wound on an insulating form around the permanent magnet. 

This signal is of course an AC signal in synchronization with the light 
chopping frequency and is used for the lock in amplifier. The coil is 
shielded with layers of Netic and Conetic metal to eliminate any 
electrical noise from being picked up in the reference circuit since 
it operates in the close vicinity of the chopper motor. 

The motor used to drive the chopper has been a problem to operate 
in the evacuated spectrometer. The 3600 rpm synchronous motors used 
originally were purchased from Bodine. The bearing grease was replaced 
with a mixture pf molybdenum disulfide and diffusion pump oil which has 
been used successfully in this laboratory in vacuum conditions. The 
motors, however, would not hold up for long periods of time. The heat 
generated by the motors is thought to be the cause of the problem. In 
a vacuum very little heat is conducted away from the motor, even though 
a large ground strap is attached to the motor mount for this purpose. 
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and it is thought that this heat led to the destruction of the motor 
directly or else vaporized the diffusion oil lubrication leading to the 
destruction of the bearings. Occasionally the fields would burn out 
from overheating before the bearings stopped the motors. Operating the 
motor at a lower voltage did not help the problem nor did other lubri- 
cants including the grease shipped in the motor. After several motors 
were tried lasting from a few days to three months depending on the 
usage, a surplus motor was obtained which was apparently not synchronous 
but very constant in rpm. This motor was tried because of its avail- 
ability and precision bearings. The motor is also totally enclosed so 
that no lubrication can be lost through vaporization. The mount was 
adapted slightly to accept the new motor which was fortunately almost 
the same physical size. Thus far this motor has operated without 
failure for six months under heavy usage. The original lubrication in 
the motor was left unchanged. It is thought that the precision bearings 
and the total enclosure of the motor allow the motor to operate properly 
even in vacuum conditions. 

3.9 Detector Pumping System 

The PbS and InSb detectors from Santa Barbara Research Corporation 
used in this experiment operate at liquid nitrogen temperatures and are 
contained in vacuum insulated dewars. A pumping system consisting of a 
mechanical forepump, a diffusion pump and a cold trap was constructed 
for evacuation of the vacuum jacket. One inch copper line was run from 
the pumping system to both dewars which are mounted in either end of 



84 


the detector section of the spectrometer, and short lengths of rubber 
hose are used to connect from the copper line to the dewars. This 
pumping system is capable of vacuums of .2 to .4 microns as measured 
on a Mercury McLeod Gauge back near the pump. The pressure is the 
dewars is probably somewhat greater. The dewars pumped with this 
system and closed off will hold liquid nitrogen for about 30 hours 
and can be refilled without repumping for 3 to 4 days. After this time, 
however, the liquid nitrogen hold time is reduced due to leakage and 
outgassing, to the point that the dewars have to be warmed up and be 
re-pumped. It is not possible to pump continuously with this system 
since the liquid nitrogen in the dewars produce a better vacuum than 
the pumps can achieve and the pumps simply backstream into the dewars. 
For this reason small one liter per second ion pumps from Ultek Corpo- 
ration were installed directly on the dewars so they could be pumped 
continuously even when they contained liquid nitrogen. The original 
pumping system is used to rough pump the dewars and ion pump manifold 
down to less than a micron and then is valved off and the ion pump is 
turned on. By pumping continuously with the ion pumps the liquid 
nitrogen hold time was increased to 50 hours. Also repumping is not 
necessary so that the dewars never have to be warmed up. The ion pump 
and manifold is shown in Figure 20 and is constructed entirely of 
stainless steel with heliarced joints to hold a high vacuum. A stain- 
less diaphragm is connected to the manifold which expands against a 
microswitch triggering an alarm if the dewar should warm up and come 
up to atmospheric pressure. 




Figure 20. Detector Dewar Pump 
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3. 10 Electronics 

Both the PbS and InSb detector signals were amplified using a 
Princeton Applied Research Type A preamplifier and HR-8 lock in 
amplifier. The amplified signal was recorded on a Texas Instrument 
chart recorder. Records were made simultaneously on chart paper and 
on magnetic tape for computer assisted analysis. In these cases the 
monitor signal from the amplifier was routed to an analog to digital 
converter and put on tape in digital form, using a PDP-11 computer. 
The software and some of the hardware for this data taking system was 
developed at The University of Tennessee Molecular Spectroscopy 
Laboratory. 



CHAPTER IV 


DATA AND ANALYSIS OF DATA 
4.1 Analysis Techniques 

Practically all of the data taken was digitally recorded on 
magnetic tape in addition to analog recordings on chart paper. For 
frequency measurements the data records included calibration lines of 
known frequency. The serial data point numbers of the known and 
unknown lines were determined from the magnetic tape record by displaying 
the data on a Tennecomp display scope using a PDP-11 computer. Software 
developed at the University of Tennessee Molecular Spectroscopy Labora- 
tory displays and moves the data across the screen face. The frequency 
coordinate can be expanded at will, and a cursor can be moved to locate 
the line center. The software keeps track of the serial data point 
number at the position of the cursor, and the data point number of the 
line center located by the cursor can be output to a teletype or line 
printer record, along with an identification of the line measured. The 
data point number of each line was measured three times, and the results 
averaged. The reproducability on the CO calibration lines was about 
one to two data points with the frequency interval between data points 
about 0.00035 cm -1 . For the H 2 and D 2 lines the reproducibility was 
somewhat poorer, due in part to a weaker signal intensity and also due 
to noisier and in some cases very weakly absorbing lines. The reproduc- 
ibility of the H 2 and D 2 measured lines varied from about 2 to 6 data 
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point numbers. The known frequencies and their serial data point 
numbers were fitted either to a quadratic or cubic equations of the 
form, 


v = A + BX + CX 2 + , 


where v is the frequency of the line and X is the data point number. 

A, B and C the coefficients were determined by a least squares fit 
Basic program using the PDP-11 computer. The program carries seven 
significant figures which was adequate after subtracting a constant 
from all the frequencies of the lines. These coefficients were used 
to determine, the unknown line frequencies. 

For line width measurements the particular line to be measured 
was recorded several times on magnetic tape. Since the deconvolution 
program as well as the response function program needs block numbers 
bracketing the data of interest (each block contains 256 data points) , 
the data from magnetic tape was plotted in a compressed form using the 
Calcomp plotting system and the UTCC, IBM 360/65 system. with the block 
numbers indicated on the plots. The block number data can then be 
quickly determined for use in the appropriate program. 

The digitally recorded data were invaluable for line width 
measurements since it eliminates the tedium and inaccuracy of transferring 
measurements from chart paper to a form which can be read by the computer. 
The speed of making frequency measurements is also increased using the 
displayed data. In both cases magnetic tape provides a more compact 
and more readily accessed data source. 
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4.2 Effect of the Electric Field on Line Position 

Before drawing any conclusion concerning the pressure shifts of 
H 2 and D 2 spectral lines it was necessary to determine the effect of 
the electric field on the line position. The lines are usually measured 
using the maximum electric field possible at a particular pressure in 
order to obtain the most intense, easily measurable line. So in addition 
to varing the pressure it is convenient to also vary the field in 
measuring the pressure shifted frequencies. To observe any electric 
field shift which might occur, the Qj (0) and Q^(l) lines of H 2 were 
recorded at a constant pressure of 19.22 Amagat at several field inten- 
sities from 160,000 to 50,000 volts/cm. The results of these measurements 
shown in Table VII and Figure 21 indicate that, within experimental error 
of about .003 cm -1 and within the range of fields used, no field dependent 
frequency shift of the H 2 lines appears. 

4.3 Pressure Shifted Frequencies of H 2 

The pressure dependent frequencies of the Q 1 (0) through Qj (3) 
spectral lines of H 2 were measured both as a check on the performance 
of the 5-meter Littrow spectrometer and because more accurate measurements 
could be made than in previous work. The pressure range used was from 
about 7 to 28 Amagats with three measurements made at each pressure 
for each line. Frequency measurements were made using the first overtone 
band of CO for calibration lines and the published frequencies of these 
lines given by Rao (15). With the 31.6 groove per millimeter grating it 
was possible to choose calibration lines either in the same or adjacent 
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TABLE VII 

FREQUENCY OF H LINES AS A FUNCTION 
OF THE ELECTRIC FIELD INTENSITY 
AT CONSTANT PRESSURE 1 


• J Ql 
Field 

(volts/cm) 

<°) , 
Frequency^ 

(cm -1 ) 

Field 

(volts/cm) 

V°> 3 

Frequency 

(cm -1 ) 

159,000 

4161.1250 

164,000 

4155.1961 

150,000 

4161.1277 

150,000 

4155.1973 

135,000 

4161.1267 

135,000 

4155.1970 

121,000 

4161.1239 

122,000 

4155.1971 

108,000 

4161.1301 

109,000 

4155.2000 

94,000 

4161.1301 

93,500 

4155.1983 

80,900 

4161.1289 

79,900 

4155.1977 

67,200 

4161.1292 

67,200 

4155.1985 

51,200 

4161.1248 

51,200 

4155.1971 

ii i 53 x : 


^Density = 19.22 Amagat 

Average Frequency = 4161.1274 ± *0024 cnT 1 
3 

0 Average Frequency = 4155.1966 ± .0037 cm -1 
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order as that of the unknown line. The 12-meter path of the long path 
cell was used to hold 5 millimeters pressure of CO and the beam switching 
mirror described in Section 3.7 was used to switch from the Stark cell 
beam to the calibration beam. Both signals were detected using the 
cooled lead-sulfide detector and recorded on chart paper magnetic tape. 
The Qj (1) lines were so close together that they could be measured 
together in the 24th order of the grating along with three calibration 
lines also in 24th order. The calibration lines were fitted to a 
quadratic equation in position using the least squares regression 
technique discussed in Section 4.1 and the unknown line frequencies 
calibrated from this quadratic equation. The Q x (2) and Q x (3) lines 
also occuring in 24th order were calibrated by linear extrapolation 
from two very close CO calibration lines occuring on either side of 
each line in 23rd order. 

Each line was measured three times at each pressure. The three 
determinations were averaged and a standard deviation calculated to 
give an indication of the reproducibility of the measurements. The 
resultant frequencies for each line were fitted to a straight line 
using a least squares method. The intercept of this line gives the 
zero pressure frequency of the line of interest and the slope gives 
the pressure shift. Tables VIII to XI show the observed and calculated 
frequenices of the lines measured, and Figures 22 to 25 show the best 
fit straight line through the data points. Table XII gives the zero 
pressure frequenices and the shifts of the H 2 lines. The slopes of the 
lines measured can be compared in Figure 25. Table XIII compares the 



TABLE VIII 

OBSERVED AND CALCULATED FREQUENCY OF THE Q (0) , H 2 LINE AS A FUNCTION OF DENSITY 


Density 

(Amagats) 

Measured Frequency 
(cm -1 ) 

Calculated Frequency 
(cm -1 ) 

Observed - 
Calculated 

.28.21 

. 4161.1014 ± .0014.: 

4161.1026 

-.0012 

25.23 

4161.1084 ± .0008 

4161.1090 

-.0006 

22.23 

4161.1153 ± .0019 

4161.1155 

-.0002 

19.22 

4161.1256 + .0019 

4161.1219 

+.0037 

16.20 

4161.1290 ± .0013 

4161.1284 

+.0006 

13.17 

4161.1343 ± .0019 

4161.1349 

-.0006 

10.12 

4161.1399 ± .0004 

4161.1415 

•V, 00.16 

7.07 

4161.1481 ± .0020 

4161.1480 

+.0001 


to 

04 


DENSITY" lAMftGflT) 



TABLE IX 

OBSERVED AND CALCULATED FREQUENCY OF THE QjCl), H 2 LINE AS A FUNCTION OF DENSITY 


Density Measured Frequency Calculated Frequency Observed- 

(Amagats) (cm* 1 )' (cm* 1 ) Calculated 


28.21 

4155.1638 

+ 

,0011 

4155.1660 

-.0022 

25.23 

4155.1750 

± 

.0009 

4155.1753 

-.0003 

22.23 

4155.1862 

+ 

.0021 

4155.1847 

+.0015 

19.22 

4155.1956 

± 

.0007 

4155.1941 

+.0015 

16.20 

4155.2051 

± 

.0009 

4155.2036 

+.0015 

13.17 

4155.2127 

± 

.0025 

4155.2131 

-.0004 

10.12 

4155.2220 

± 

.0010 

4155.2226 

-.0006 

7.07 

4155.2311 

± 

.0011 

4155.2322 

-.0011 




FREQUENCY (CM-*) 

Figure 23.- Q^l), H 2 Line Position versus Density 



TABLE X 


OBSERVED AND CALCULATED FREQUENCY OF THE Q x (2) , H 2 LINE AS A FUNCTION OF DENSITY 


Density 

(Amagats) 

Measured Frequency 
(cm -1 ) 

Calculated Frequency 
(cm - 1 ) 

Observed- 

Calculated 

28.21 

4143.4139 ± .0021 

4143.4136 

-.0007 

25,23 

4143.4205 ± .0020 

4143.4207 

-.0002 

22.23 

4143.4281 ± .0030 

4143.4269 

+.0012 

19.22 

4143.4344 ± .0009 

4143.4331 

+.0013 

16.20 

4143.4384 ± .0023 

4143.4393 

t . 0009 

13.17 

4143.4447 ± .0048 

4143.4455 

-.0008 

10.12 

4143.4521 ± .0037 

4143.4518 

-.0003 
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FREQUENCY (CM -i ) 


Figure 24. Q/2), H 2 Line Position versus Density 



TABLE XI 


OBSERVED AND CALCULATED FREQUENCY OF THE Q x (3) , H 2 LINE AS A FUNCTION OF DENSITY 


Density 

(Amagats) 

Measured Frequency 
(cm -1 ) 

Calculated Frequency 
(cm" 1 ) 

Observed- 

Calculated 

28.21 

4125.8326 + .0016 

4125.8332 

- . 0006 

25.23 

4125.8410 ± .0025 

4125.8386 

+.0024 

22.23 

4125.8441 ± .0076 

4125.8440 

+.0001 

19.22 

4125.8473 ± .0034 

4125.8495 

-.0022 

16.20 

4125.8518 ± .0032 

4125.8550 

-.0032 

13.17 

4125.8638 ± .0016 

4125.8605 

+.0033 


U3 

KD 


i 
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TABLE XII 

CALCULATED ZERO PRESSURE FREQUENCY AND 
PRESSURE SHIFT OF H 2 LINES 


Line 

Frequency 
(cm - *) 

Pressure Shift 
(cm _1 /Amagat) 

Qj CO) 

4161.1632 ± .0017 

-.00215 ± .89xl0 -4 

QjCl) 

4155.2543 ± .0014 

-.00313 ± ,75xl0" 4 

Qj(2) 

4143.4725 ± .0013 

-.00205 ± . 64xl0 -4 

Ql(3) 

4125.8845 ± .0048 

-.00182 ± . 23x10" 3 
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TABLE XIII 


h 2 frequencies and pressure shifts measured in this experiment compared to others 


Line 

Frequency Frequency^ 1 

Measured this Measured by 

Experiment (cm -1 ) Brannon (cm - 

Frequency 2 
Measured by 
Stoicheff (cm - 

Frequency 3 
Measured by 
Fink et al. (cm - 

31 ■■ ' ,„Ji r 

Frequency ^ 
Measured by 

*) Foltz et al. (cm -1 ) 

tyo) 

4161.1632 

4161.158 

4161.134 



(yi) 

4155.2543 

4155.240 

4155.202 

4155.2575 

4155.259 

Qj(2) 

4143.4725 

4143.449 

4143.387 

4143.4688 


V 3) 

4125.8845 

4125.867 

4143.832 

4125.8718 



Pressure Shift Pressure Shift 

Line Measured this Measured by 

Experiment Brannon 

(cm -1 /Amagat) (cm -1 /Amagat) 

Pressure Shift 
Measured by 
Fink et al. 
(cm - VAmagat) 

Pressure Shift 
Measured by 
Foltz et al . 
(cm -1 /Amagat) 


Q x (0) -.00215 


-.0018 


- . 00234 


Qj(l) -.00313 


-.0029 

-.00240 

-.00318 


Ql (2) -.00205 


-.0015 


-.00202 


Qj (3) -.00182 


-.0018 


-.00185 


Electric Field Induced Technique 
2 

Raman Technique 

3 

Quadrupole Technique 

4 

Electric Field Induced and Raman Technique 
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frequencies and shifts measured in this experiment with those measured 
by Brannon (5), Stoicheff (16), Fink et al. (17) and Foltz et al. "(18). 

The frequencies measured in this experiment are shifted with 
respect to those quoted by Brannon and Stoicheff, but agree very well 
with the values measured by Fink et al. and Foltz et al. The calibration 
techniques used and the resolution available in this experiment were 
better than that available to either Brannon or Stoicheff. Their 
quoted precisions were ,02cm" 1 and. 01cm" 1 respectively whereas this 
experiment is able to quote a precision of about .002cm" 1 . The 
agreement of the linear shift coefficients measured. in this experiment, 
with those measured by Foltz et al. is very good. 

Both the relative and absolute accuracy of the frequency 
measurements made in this experiment are felt to be very good. The 
direct calibration techniques used and the resolution available were 
excellent and make highly accurate measurements possible, 

4.4 Pressure Shifted Frequencies of P 2 

The pressure shifted frequencies of the Q 1 (0) through Q(4), 
lines were measured since more precise measurements were possible than 
in previous work. The methods used were essentially the same as 
employed for H 2 described in the previous section. All of the lines, 
occuring in 17th order of the 31.6 groove per millimeter grating, were 
run together with 8 calibration lines of CO occuring in 24th order, 
interspersed among the Q-branch lines. A cubic equation was fitted to 
the calibration lines using a least squares technique and the unknown 
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lines calibrated from this curve. Again as for H 2 a straight line was 
fitted to the line positions as a function of density to obtain a zero 
pressure intercept frequency and pressure shift coefficient. Tables 
XIV through XIX : show the measured pressure shift data and the calcu- 
lated values from the best fitted straight line. Figures 27 through 31 
show the best fit line drawn through the measured data points and 
Figure 32 compares the slopes of the D 2 Q-branch lines. 

A comparison of the frequencies and pressure shifts measured in 
this experiment as compared to measurements made by Brannon (5) and 
Stoicheff (16) is shown in Table XX. Again the frequencies measured 
in this experiment are shifted from those measured by Brannon or Stoicheff 
but agree more closely with Brannon. It is felt that the frequencies 
measured here are the more accurate of the two. The accuracy of the 
frequencies and shifts measured in this experiment was demonstrated 
by the H 2 data discussed in the previous section. Again superior 
resolution and the calibration techniques used are felt to contribute 
to more accurate measurements 

4 . 5 Field Induced Branch of N 

The unresolved Q branch of N 2 shown in Figure 33 was observed but 
was too weak to make any new quantitative measurements . The peak 
absorption observed for the Q branch was 3.5% with no zero off set at 
an electric field intensity of 220,000 volts/cm. This measurement was 
made with no polarizer in the beam. The closely spaced lines of the Q 
branch were not expected to be resolved at the 28.21 Amagat pressure 



TABLE XIV 


OBSERVED AND CALCULATED FREQUENCY OF THE Q x (0), D 2 LINE AS A FUNCTION OF DENSITY 


Density Measured Frequency Calculated; Frequency Observed- 

(Amagat) (cm -1 ) (cnf 1 ) Calculated 


28.21 

2993.5724 

± 

.0013 

2993.5686 

+.0038 

25.23 

2993.5735 

± 

.0029 

2993.5743 

-.0008 

22.23 

2993.5773 

+ 

.0030 

2993.5801 

-.0028 

19.22 

2993.5825 

± 

.0002 

2993.5859 

-.0034 

16.20 

2993.5912 

± 

.0042 

2993.5917 

-.0005 

13.17 

2993.6014 

+ 

.0029 

2993.5976 

+.0038 
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FREQUENCY (CM _1 ) 


Figure 27. Q 1 (0) ; D 2 Line Position versus Density 
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TABLE XV 


OBSERVED AND CALCULATED FREQUENCY OF THE Q : (1) , D 2 LINE AS A FUNCTION OF DENSITY 


Density Measured Frequency Calculated Frequency Observed- 

(Amagat) (cm -1 ) - (cm -1 ) Calculated 


28.21 

2991.4605 ± 

.0009 

2991.4592 

+.0013 

25.23 

2991.4652 ± 

.0009 

2991.4646 

+.0006 

22.23 

2991.4673 ± 

.0008 

2991.4699 

-.0026 

19.22 

2991.4747 ± 

.0010 

2991.4754 

-.0007 

16.20 

2991.4809 ± 

.0005 

2991.4808 

-.0001 

13.17 

2991.4880 ± 

.0021 

2991.4863 

+.0017 

10.12 

2991.4914 ± 

.0020 

2991.4918 

-.0004 



DENSITY (RMflGHT) 



Figure 28. Q^l), d 2 Line Position versus Density 
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TABLE XVI 


OBSERVED AND CALCULATED FREQUENCY OF THE Qj (2) , D 2 LINE AS A FUNCTION OF DENSITY 


Density 

(Amagat) 

Measured Frequency 
(cm" 1 ) 

Calculated Frequency 

Observed- 

Calculated 

28.21 

2987.2406 ± .0035 

2987.2395 

+.0011 

25.23 

2987.2454 ± .0018 

2987.2458 

- . 0004 

22.23 

2987.2503 ± .0016 

2987.2521 

-.0018 

19.22 

2987.2584 ± .0005 

2987.2585 

-.0001 

16.20 

2987.2642 ± .0018 

2987.2648 

-.0006 

13.17 

2987.2737 ± .0023 

2987.2712 

+.0025 

10.12 

2987.2786 ± .0019 

2987.2776 

+.0010 

7.07 

2987.2825 ± .0017 

2987.2841 

-.0016 





TABLE XVII 


OBSERVED AND CALCULATED FREQUENCY OF THE Qj (3) , D 2 LINE AS A FUNCTION OF DENSITY 


Density 

(Amagat) 

Measured Frequency 
(cm" 1 ) 

Calculated Frequency 
(era -1 ) 

■ Observed- 
Calculated 

28.21 

2980.9465 ± .0032 

2980.9450 

+ .0015 

25.23 

2980.9496 ± .0068 

2980.9501 

-.0005 

22.23 

2980.9537 ± .0045 

2980.9552 

-.0015 

19.22 

2980.9603 ± .0047 

2980.9603 

..0000 

16.20 

2980.9642 ± .0045 

2980.9655 

-.0013 

13.17 

2980.9724 ± .0112 

2980.9707 

+.0017 
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TABLE XVIII 


OBSERVED AND CALCULATED FREQUENCY OF THE 0^4), LINE AS A FUNCTION OF DENSITY 


Density Measured Frequency Calculated Frequency Observed- 

(Amagat) (cm -1 ) (cm -1 ) Calculated 


28.21 

2972.5710 

± 

.0086 

2972.5726 

-.0016 

25.23 

2972.5799 

± 

.0038 

2972.5779 

+.0020 

22.23 

2972.5838 

± 

.0072 

2972.5832 

+.0006 

19.22 

2972.5884 

± 

.0108 

2972.5885 

-.0001 

16.20 

2972,5925 

+ 

.0051 

2972.5938 

-.0013 

13.17 

2972.5996 

± 

.0027 

2972.5992 

+.0004 
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DENSITY (RMRGRT) 



Figure 31, Q. 1 (4) J D 2 Line Position versus Density 
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TABLE XIX 

CALCULATED ZERO PRESSURE FREQUENCY AND 
PRESSURE SHIFT OF D 2 LINES 


Line Frequency Pressure Shift 

(cm -1 ) (cm -1 /Amagat) 


QiCO) 

2993.6230 

± .0059 

-.00193 

± 

.28x10" 3 

Q x (!) 

2991. S099 

± .0020 

-.00180 

± 

.10x10" 3 

Q x (2) 

2987.2990 

± .0015 

-.00211 

± 

.80xl0 -4 

Qj(3) 

2980,9932 

± .0026 

-.00171 

± 

.12x10" 3 

Q z (4) 

2972.6224 

± .0026 

-.00176 

± 

. 12x10" 3 
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TABLE XX 

d 2 FREQUENICES AND PRESSURE SHIFTS MEASURED 
IN THIS EXPERIMENT COMPARED TO OTHERS 


Line 

Frequency 
Measured this 
Experiment (cm -1 ) 

Frequency 
Measured by 
Brannon (cm - *) 

Frequency 
Measured by 
Stoicheff (cm 

Q 1 (0) 

2993.6230 

2993.600 

2993.548 

Q x CO) 

2991.5099 

2991.486 

2991.446 

Ql(2) 

2987.2990 

2987.268 

2987.230 

QiC3) 

2980.9932 

2980.967 

2980.877 

Ql(4) 

2972.6224 

2972.597 

2972.557 


Line 

Pressure Shift 
Measured this 
Experiment 
(cm“ VAmagat) 

Pressure Shift 
Measured by 
Brannon 
(cm“ VAmagat) 

Qi (0) 

-.00193 

-.0017 

Qi(l) 

-.00180 

-.0016 

QiC2) 

-.00211 

-.0017 

Ql(3) 

-.00171 

-.0014 

Qi(4) 

-.00176 

-.0020 
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used to observe the lines; however the wider spaced S and 0 branch lines 
could be used for quantitative measurements. These lines are however 
four to five times less intense than the Q branch lines and could not 
be detected. 

Some problems were encountered in trying to detect the N 2 lines . 

The N 2 frequencies lie in a region where C0 2 is strongly absorbing. 

The flushing cover described in Section 3.6 was successful in removing 
practically all of the C0 2 in the light path from the source to the 
Stark cell, however as observed in Figure 33 on page 119, some small 
amount remains. The bottom trace show the C0 2 background in the light 
path with the electric field off. The zero has been offset to increase 
the recorded signal. The top trace has been offset slightly less, and 
shows the N 2 , Q branch after the field has been turned on and the same 
region rescanned. The two traces are aligned only approximately. Most 
of the C0 2 seen here is in the flushing cover. These lines were used 
to calibrate the frequency of the N 2 , Q branch. A greater problem 
however was encountered with C0 2 absorption in the Stark cell itself. 

It was found that C0 2 appeared slowly and grew within the cell over a 
period of 24 to 36 hours. Repeated evacuation and refilling did not 
seem to slow down the reappearance of the C0 2> Since the cell: is 
pressurized to 4S0 p.s.i. the C0 2 appeared as a broad absorption envelope, 
and until the problem was discovered, masked the small Q branch absorption. 
The source of the C0 2 is not known but it is possible that it could be 
generated somehow by the plexiglas or other materials used in the cell. 
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The center of the absorption of the Q branch was measured to be 
2328.8 cm -1 . From molecular constants given by Herzberg (3) the 
frequencies and relative intensities of the N 9 lines were calculated 
and show the center of the Q branch should occur at a frequency of 
2329.67 cm -1 corresponding to the (7) line. This is a total shift 
of about -0.9 cm' 1 or about -.03 cm -1 /amagat. A measurement of the 
full width at half maximum of the observed absorption revealed a width 
of 1.0 cm -1 indicating that the (^(3) through Q 1 [9) lines were contri- 
buting significantly to the absorption. 

4 . 6 Polarizability Constants of H 2 

The polarizability of H 2 was determined using the linear absorption 
approximation given in equation (55) . The integrated intensity of the 
(0) line, which depends only on the polarizability constant a, was 
measured at several pressures and electric field intensities with a 
6-plate glass polarizer in the beam to allow only the parallel polarized 
component of the incident radiation through the cell. No instrument 
function correction is needed for this measurement since according to 
Brannon (5) the integrated area is approximately independent of the 
instrument function. The area is then normalized to the value of the 
100 percent absorbing point which was determined by pressurizing the 
Stark cell with about 125 p.s.i. of carbon monoxide in order to make 
the weak overtone line in that region 100 percent absorbing. This value 
is adjusted using the efficiency of the polarizer determined from the 
maximum absorption of Q 2 (1) line at 31.22 amagat and a field of 193,000 
volts/cm. In this particular case the area was measured using a polar 
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planimeter directly from the recorded data on chart paper. The 
normalized area is then multiplied by the wave number interval per inch 
determined from calibration lines to give the resultant integrated 
absorption in units of cm -1 . A total of 20 individually measured areas 
were used to determine a. The large number of measurements were used 
because of the inherent large error in measurements of the area. The 
individual measurements were scattered randomly about an average value. 
Since y makes only about a 3 percent or less, contribution to the 
intensity of the Q branch lines it was not possible to determine this 
constant from the Q branch lines. The S^(l) line was used to determine 
the polarizability y for H^. The strength factors for the S-branch 
lines depend only on the polarizability constant y, and S^l) line is 
the strongest of the S lines. The line is never the less still very 
weak and so the line was run 15 times at 28.21 amagat and the maximum 
field possible, to get a measurable absorption. The area was normalized 
as before and converted to an integrated absorption in units of cm" 1 . 

Equation (57) and (55) along with the appropriate Boltzmann factor 
from Table V on page 38 and the appropriate strength factor from Table IV 
on page 32 were used to calculate the constants ct and y given in Table 
XXI. A standard deviation of the individual measurements for both 
determinations is also given in Table XXI, A comparison of the polariz- 
ability constants measured by Brannon (5) and Foltz (15) shows good 
agreement for the constant y between all experiments. The agreement for 
a is not as good. To check the value of a calculated in this experiment. 
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TABLE XXI 


POLARIZABILITY 

CONSTANTS OF H 2 

a Measured 
this Experiment 
(cm 3 ) 

Y Measured 
this Experiment 
(cm 3 ) 

0.802 x 10" 25 ± .049 x 10“ 2S 
0.805 x 10" 251 

0.772 x 10 -25 ± .035 x 10~ 25 

a Measured 
by Brannon 
(cm 3 ) 

Y Measured 
by Brannon 
(cm 3 ) 

0.97 x 10‘ 25 ± 3% 

0.81 x 10" 25 ± 5% 

a Measured 
by Foltz 
(cm 3 ) 

Y Measured 
by Foltz 
(cm 3 ) 

1.04 x 10" 25 

0.79 x 10" 25 

^Measured from the deconvoluted 

data of this experiment. 
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the integrated absorption data taken from the absorption coefficient 
curves of the {^convoluted Qj (P) , H 2 line were used to. recalculate ;'a. 

The: use; of this data removes -then restrict ion of the linear absorption- ... 
approximation: As seen inofable XXI the recalculated, value of- a usings 
the deconvplutedydata is virtually they same^as - the previously^ calculated 
Value. The reason for the discrepancy is the values -calculated for the 
different experiments is not known. 

4.7 Deconvolution of Spectral Lines 

In order to determine more accurately the true width of a spectral 
line it is necessary to subtract out the influence of the spectrometer 
on the observed line profile. This influence is mainly due to finite 
width slits but also includes such effects as optical abberations and 
electronic distortion. If we can determine accurately enough the 
influence of the spectrometer which we will call the response function, 
then it can be removed by solving the convolution integral 


Fo(v) 



-00 


R(v-xOF T (\Odv' 


> 


where Fo is the observed fractional absorption, R .is the response 
function and F^, is the true fractional absorption. Determination of 
the response function is discussed in detail in Section 4.8. To solve 
for the response function using computer techniques the convolution 
integral is usually written as 


(Fo)i - Z R u (P T ) t ■ 



The commonly used solution to this equation is the Van Cittert 
approximation. This is an interative process and the K+l approximation 
is according to Jansson et al. (19). 


, (K+l) _ (K) 

(F t ). - (F T ) i 


n 


(Fo) - l 
£=1 




Jansson has modified this solution to assure convergence and good 
physical results by introducing a variable parameter k, multiplying 
the applied correction factor to control the amount of correction used 
in each interaction. This prevents over correction and rapidly 
increasing oscillations in the solution which result from over correction. 
Jansson used a parameter which also was varied with the value of the 
fractional absorption at the point being corrected. The parameter was 


K. 


1 



(Vi 


(K) 



9 


where x is a constant with a value between zero and one. This factor 
o 

is in effect a weighting factor which varies the constant correction 
factor parameter from a maximum at 50 percent absorption to zero at 
both 0 and 100 percent absorption. This factor was not used in the 
deconvolution programs used in this work since it seemed possible that 
the line shape of the final result could be modified by the inclusion 
of the factor. 

Jansson' s final expression for the K+l approximation (substituting 
the constant k for the weighting factor V is 
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cvr 1 ■ CF I )f ) * K: 


< F °L -Jj m« (I W 


(K+l) 


n 

■ I 

JL-i 


l (R) u CF T )f>f , 


where it can be noted that the calculated (F^)^ values are used 
immediately as they are obtained. The deconvolution program works 
in four loops. Loop zero uses a polynomial fit to smooth the input 
data a variable number of times without deconvolution. Loop one 
deconvolutes and then smoothes the input data a variable number of 
times with an input smoothing parameter k. Loop two repeats the same 
process with a different jc if desired. Loop three deconvolutes without 
smoothing, a variable number of times again with a different k if 
desired. Jannson’s program was modified to deconvolute only every 
seventh point in loop one and every fourth point 1 in loop two with the 
points in between obtained by interpolating, thus increasing the speed 
of the program. In the H 2 data used, one initial smooth was used with 
three, four and five passes used respectively in loop one, two and three. 
The k values used were 0.3, 0.5 and 1.0 respectively in each of the 
three loops. The complete program for the deconvolution process is 
listed in Appendix B and C, and includes Jansson’s modified program 
in the subroutines. The widths of the lines were observed to vary some 
what due to noise, which made several determinations* of each line width 
necessary. In order not to have to repeat the time consuming deconvolu- 
tion process for each individual measurement, the program was written 



127 


to be able to take several measurements of the same line at the same 
conditions and align and average them so that the deconvolution is 
performed only once on the averaged data. This has the additional 
advantage of reducing the noise on the lines which improves the perfor- 
mance and results of the deconvolution process. The individual 
measurements are smoothed once and the base line subtracted out. The 
100 percent absorbing value and polarizer efficiency is used to 
convert the data to a fractional absorption and the center of each 
line measurement is found at the value of the 50 percent peak absorbing 
value. The progress of the program is checked at various points by 
plotting the data using the Calcomp plotter and the UTCC IBM 360/65 
system. The deconvoluted data are converted to an absorption coefficient 
curve and plotted in the final steps. The program measures the half 
width of the absorption coefficient curve and also integrates the area 
under the curve. Figures 34 and 35 show the results of the averaging 
process for the (0) line of H 2 , before and after deconvolution. 

Figure 36 shows the absorption coefficient curve calculated from the 
deconvoluted line. 

4.8 Determination of the Response Function 

In order to determine the response function used in the 
deconvolution program, a line was chosen in the very close vicinity 
of the lines of interest. The line chosen was negligibly broadened 
by pressure effects, so that the line width is due solely to the Doppler 
broadening> instrumental broadening and to a negligible degree to natural 
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broadening, pressure broadening was eliminated by admitting CO into 
the spectrometer tank. With a 72 meter path, 0.6 millimeters of carbon 
monoxide was sufficient pressure to use the P(21) line in the first 
overtone with no observable pressure broadening. This line was chosen 
for the deconvolution of the Q^(0) and Q^(l) lines since its 
frequency was between the two lines of interest and was in the same 
order as the H 2 lines. The CO line was observed at a scan rate and 
slit width corresponding to the data to which it would be applied. The 
line was observed several times and recorded on magnetic tape as well 
as on chart paper. The portions of the runs containing the lines were 
then smoothed and plotted from the magnetic tape record using the IBM 
Calcomp plotting system. Figures 37 and 38 show the raw data and the 
smoothed data plotted out by the Calcomp plotter. The same process 
of averaging several measurements was used for the response function 
as was used for the data in the deconvolution process. After the base 
line was subtracted and the data converted to fractional absorptions, 
the line centers were used to align the several measurements and the 
average was found. The result of this process is shown in Figures 39 
and 40. This averaging process results in a significant reduction of 
the noise on the response function. The observed line profile is a 
convolution of the actual response function and the natural Doppler 
width. The Doppler profile is gaussian and can be calculated from the 
formula, 



A(v) = 
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where is the frequency of the observed line, m is the mass of the 
CO molecule, K is Boltzmann's constant and c is the speed of light. 

The deconvolution program may now be used to deconvolute the observed 
line profile, removing the calculated Doppler contribution and leaving 
the response function which describes the instrumental broadening. 

These steps are shown in Figures 41 and 42 where the scale of the 
response function has been increased to observe more readily the effects, 
of the deconvolution. The excess base line is discarded and then the 
area under the line profile is computed. Each point was then divided 
by the area under the response function curve to normalize the area to 
unity, and the results are punched on cards for use in the deconvolution 
program. The final response function used for more of the work here 
contained 352 points. The final response function is compared to a 
calculated Gaussian of the same half width and height in Figure 43. 

It is observed that the result is assymetric when compared to the 
Gaussian. The cause of this assymetry is not known but may possibly 
be due to abberations in the spectrometer optics or the use of straight 
entrance and exit slits. Curved slits are being built to replace the 
straight slits and the effect of the replacement of the slits on the 
symmetry of spectral lines should be interesting. Appendix B lists 
the computer program used for the determination of the response function 
and gives further specific details of the program. 

4.9 Width of the H 2 , Q 1 (0) and Q t (l) Lines 

The widths of the H 2 , Q branch lines have been measured and 
calculated by Hunt et al. (20). The paper gives references to the 








Figure 43. Response Function Compared to a Gaussian 
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theory of Van Kranendonk (21) and other later expositions pertinent to 
the H 2 line widths. Brannon (5) observed a possible dependence of the 
line width on the electric field intensity. Hunt et al. observed no 
field dependence on the width but offered no data on that effect. For 
this reason two of the stronger Q branch lines, the Qj (0) and Q^fl) were 
observed at a constant density of 28.21 amagat and varying field 
intensities. Table XXII and Figure 44 show the data taken for this 
experiment and show no field dependence of the width of the line to 
within experimental error. The half widths are measured from the 
absorption coefficient curves produced by the deconvolution program 
described in Section 4.7. 

The absorption coefficient curve is approximately Lorentzian in 
shape as it should be for pressure broadened lines. This is observed 
in Figure 45 where the Q 1 (0) , H 2 line is compared to a Lorentz function 
of the same height and half width. The broadening coefficients as a 
function of density of these two lines were measured and compared to 
those measured by Hunt et al. (20). Table XXIII and Figures 16 and 17 
show the width of the Q^O) and Q 2 (1) line as a function of density. 

The broadening coefficients measured in this experiment agree to within 
experimental error to those observed previously despite the relatively 
fewer points used in this determination. 

4.10 The Response Function at Wider Slits 

Figure 48 shows a comparison of two response functions, one obtained 
at the slit width used for the deconvolution of the lines discussed in 
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TABLE XXII 

WIDTH OF H LINES AS A FUNCTION 
OF THE ELECTRIC FIELD INTENSITY 
AT CONSTANT PRESSURE 1 


Field 

(volts/cm) 

Ql (!) 2 

Width 

(cm - l ) 

Field 

(volts/cm) 

Ql(°) 3 

' Width -5 
(cm -1 ) 

179,000 

,068. 

179,000 

.085 

160,000 

.079 

160,000 

.083 

141,000 

.060 

141,000 

.076 

122,000 

.073 

122,000 

.089 

103,000 

.065 



84,000 

.059 



65,800 

.078 



47,000 

.076 




^Density = 28.21 Amagat 
^ Average Width = .070 ± .008 cm -1 
Average Width - *083 ± ,005 cm 1 
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TABLE XXIII 


WIDTH OF H 2 LINES AS A FUNCTION OF DENSITY 


Density 

(Amagat) 

Qx CO) 

Width 

(cm -1 ) 

Density 

(Amagat) 

Qi (0) 

Width 
(cm" 1 ) 

28.21 

.083 

28.21 

.070 

25.23 

.087 

23.73 

.051 

22.23 

.075 

19.22 

.045 

19.22 

.061 

14.69 

.046 



10.12 

.041 


Broadening Coefficient 
Measured this: Experiment 
(cm -1 /Amagat) 

Broadening Coefficient 
Measured, this Experiment 
(cm" ^/Amagat) 


.0026 ± .0010 

.0014 ± .0005 


Broadening Coefficient 
Measured by Hunt et al. 
(cm"i/Amagat) 

Broadening Coefficient 
Measured by Hunt et al. 
(cm -1 /Amagat) 



.00278 


.00183 








148 


Section 4.8 and the other obtained using wider slits. The full half 
width of the narrow response function is .041 cm -1 . The wider response 
function obtained at the wider slit setting has a full half width of 
.060 cm -1 and has been rescaled to the same half width as that of the 
narrow response function and aligned at the maximum absorption point. 
This was done to compare the shapes of the two functions. It is readily 
observed that the response function observed at wider slits is still 
assymetric in the same manner as before. However the wing on the low 
frequency side of the line has dropped somewhat and appears to be more 
gauss i an. 

The response function taken at the wider slit width was used to 
deconvolute the (^(1), H 2 line observed at the same wider slit width. 

The width of the absorption coefficient curve measured at the wide 
slit width was .065 cm -1 as compared to .059 cm -1 at the narrower slit 
width. The electric field and pressure were the same for the two 
observations. The measured area at the smaller slit width was 7 percent 
larger than that measured at the wider slit setting. The total number 
of deconvolution iterations was increased from 12, which was used rou- 
tinely in most of this work using the narrow response function, to 18 
iterations using the wider response function. The result using 12 total 
iterations was essentially the same as for 18. 

The difference between the two measured widths is within the range 
of variations which was observed on individual measurements of the line 
widths. So within the experimental error encountered, the deconvolution 
program used here yields the same half widths at both slit widths used. 
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Further tests should be made to make more conclusive statements. If 
the deconvolution method used here is good enough to allow the sacrifice 
of some resolution for an improvement of signal to noise ratio then the 
time and effort expended in the acquisition of data could be reduced 
and the final results improved by the use of wider slits. The reproduc- 
ibility of line widths and areas is very much affected by the amount of 
noise in the data, and the deconvolution process itself works better 
when less noise is present in the data. An improvement in the signal 
to noise ratio would be extremely important to work similar to that 
done here where relatively weak signals are observed at high resolution. 

4,11 Molecular Constants of and 

The number of zero pressure frequencies of the Q branch lines 
determined in this experiment is not enough to determine all of the 
molecular constants of H 2 and D,,. However the differences in the 
ground and first excited states could be determined, i.e. Bj-Bq, Dq-Dj 
and Hj-Hqv Because the Q branch lines are so close together the term 
in the difference of H was not large enough to contribute, significantly 
to the equation describing the Q-branch frequencies. The Q branch lines 
were therefore fitted to the equation (see Section 2.1). 

v = v 0 ^ + (B 1 -B 0 )J(J+1) + (Dg-D^J^J+l) 2 . 

The differences in the B and D values and the intercept frequency 
were determined using a least squares method. Table XXIV shows the 
calculated constants for both H£ and ^ 2 * and the calculated values using 
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TABLE XXIV 
H 2 and D 2 .. CONSTANTS' 


V 0+l 

4161.1630 95% 

confidence 

interval 


b k b o 

-2.9570 + 

.0025 


V d i 

1.43 x 10" 3 ± 

.20 x 

10~ 3 

Observed 

Calculated 


Obs . -Cal . 

4161.1632 

4161.1630 


+.0002 

4155.2543 

4155.2547 


-.0004 

4143.4725 

4143.4723 


+.0002 

4125.8845 

4125.8845 


.0000 


D 2 Constants 



v 0+l 

2993.6223 95% 

confidence 

interval 


B r B 0 

-1.0559 + 

.0017 


D o- D i:: 

+2.95 x 10 -4 ± 

.83 x 

10“ 4 

Observed 

Calculated 


Obs. -Cal. 

2993.6230 

2993.6223 


+.0007 

2991.5100 

2991.5116 


-.0016 

2987.2990 

2987.2974 


-.0016 

2980.9932 

2980.9939 


-.0007 

2972.6224 

2972.6222 


+.0002 
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these coefficients. It should be noted that the residuals are smaller 
than the experimental error. 

The constants determined from the quadratic fit are compared 
with the appropriate data determined by Brannon (5) , Stoicheff (14) 
and Foltz et al. (18) in Table XXV. The difference in the B values 
measured here compare favorably with the previous measurements. The 
difference in the D values measured here is smaller than previous 
measurements. These previous measurements used a greater number of 
lines which would possibly increase the value of the D constants to 
counteract the inclusion of constants involving H. 

4.12 Integrated Areas of the H ? , Q 1 (0) and Q, (1) Lines 

The absorption lines observed in this experiment are pressure 
broadened and should therefore be Lorentizian in profile. The equation 
describing the absorption coefficient is 


K 


v 


S Av 

7T Av 2 +(V-Vq) 2 


(59) 


where Av is half of the full line width at half maximum absorption. v Q 
is the central frequency of the line and S is the integrated absorption 
coefficient defined as 


S = 


K d 

v v 


or in other words the area under the absorption coefficient curve. 
According to these two equations we see that given the half width and 
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TABLE XXV 

H 2 and D 2 CONSTANTS . COMPARED TO OTHER MEASUREMENTS 



; v 0+l 

B r B o 

D 0' D 1 

This work 

4161.1630 

-2.9570 

+1.43 x 10" 3 

Brannon 

4161.158 

-2.9627 

+1.82 x 10“ 3 

Stoicheff 

4161.134 

-2.9614 

+1.64 x 10' 3 

Foltz et al. 

4161.1815 

-2.9650 

+2.03 x 10- 3 






V 0-+l 

B r B o 

D 0' D 1 

This work 

2993.6223 

-1.0559 

+2.95 x 10~ 4 

Brannon 

2993.606 

-1.0575 

+4.38 x 10' 4 

Stoicheff 

2993.548 

-1.0623 

+5.9 x 10“ 4 
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and integrated area of the absorption coefficient curve we should be 
able to calculate the absorption coefficient profile. 

The Q 1 (0) and Q L (1) lines were deconvoluted in this experiment. 

After the deconvoluted profiles were converted to absorption coefficient 
curves the half widths and integrated areas were measured. The half 
widths of these lines were not observed to be a function of the electric 
field intensity, however the integrated area should, according to 
equation (57), increase as the square of the electric field intensity. 
Figure 49 clearly shows this to be the case, where the integrated areas 
of Q : (0) and Q L (1) have been measured as a function of electric field 
intensity squared at a constant density of 28.21 amagat and fitted with 
a straight line calculated using a leasts squares method. The integrated 
area should also, according to equation (57) increase linearly with the 
pressure, so that the product of the pressure and electric field 
intensity squared versus the integrated area should yield a straight 
line. Figure 50 demonstrates the validity of this argument. 

It should be emphasized that the area measurements were made' on the 
absorption coefficient curve. This is necessary since nonlinear 
absorption apparently affects the integrated areas as well as the more 
obvious effect on the half widths. A measurement of the integrated 
areas of the strongly absorbing (1) line before deconvolution versus 
the product DE 2 did not yield a straight line. 

If we assume that the integrated area S is given by the expression 

S = m DE 2 , 
s 




DENSITYwFIELD SQUARED (MOM VO 



Figure 50. Integrated Area versus DE 
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where D is the gas density, E is the electric field intensity and :i 
is the slope of the best fit line shown in Figure 50, then the value 
of S at an arbitrary density and electric field intensity can be 
computed for use in equation (59) . If we further assume that the half 
widths Av are given by the expression 

Av = m. D , 

Av 

where m is the broadening coefficient then the half width at an 
arbitrary pressure can be calculated for use in equation (59) . These 
values can then be used to calculate the absorption coefficient curve 
at the desired value of D and E. Table XXVI gives the values of the 
slopes m s and m Ay used for this purpose * The value of m g was obtained 
from the data shown in Figure 50. The value of m^ y used is the value 
measured by Hunt et al. (17) since the present work used only a limited 
amount of data to observe the broadening coefficient. Table XXVI also 
compares the calculated absorption coefficient to an experimentally 
measured coefficient. The values are actual coefficients and are not 
normalized to the path length. The calculated values compare satisfac- 
torily with the observed values. Therefore for the two lines observed 
here (Q 1 (1) and (^(0)) we can calculate the absorption coefficient profile 
at any desired electric field intensity and gas density. 
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TABLE XXVI 

CALCULATION OF ABSORPTION COEFFICIENTS 


m 


cm 


-l 


s (amagat) [VoltsT^ 


cm 



cm -1 /amagat 


(^(0) .1.32 x Id- 15 

2.78 x 10" 3 

Q l ( 1) 7.08 x 10“ 15 

1.38 x 10" 3 


K Calculated 

K Measured 

V 

V 

2.333 

2.026 

1.998 

1.831 

1.394 

1.481 

.927 

1.151 

.631 

.857 

.447 

.624 

.330 

.455 

.252 

.335 

.197 

.264 

.159 

.224 

.130 

.184 

.109 

.145 

.097 

.110 

.083 

.080 

.072 

.069 

.063 

.060 

.055 

.048 

.048 

■ .042 


^Measured by Hunt et al. (17) 


CHAPTER V 


SUMMARY AND SUGGESTIONS FOR FUTURE WORK 


5.1 Summary 

To summarize the experimental section of this work the following 
things have been achieved. 

1. The Qj (0) through Q 1 (3) , H 2 fundamental absorption line 
frequenices have been measured as a function of density. The shift 
with pressure in the 7 to 28 amagat range was found to be linear 
and zero pressure frequencies and shift coefficients are calculated 
and compared to previous values. The relative and absolute accuracy 
of the values measured are thought to be better than previous measure- 
ments because of the instrumental resolution available and the calibra- 
tion techniques used. 

2. A check on the effect of the electric field on the line 
position of the (0) and Q L (1), H 2 lines showed no shift effects 
within an experimental error of about .002 cm -1 . The field intensities 
used ranged from about 50,000 to 160,000 volts/cm. 

3. The Qj^ (0) through Q 1 (4) , D 2 fundamental absorption line 
frequencies were measured as a function of density using the same method 
as for the H 2 . The calculated zero pressure frequencies and pressure 
shifts are thought to be better than the relatively few previous measure- 
ments on these lines. 
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4. The differences in the ground and first excited state molecular 
constants B and D were found and used to calculate the zero pressure 
frequencies to less than experimental error. A comparison of these 
constants to other previously measured . values was satisfactory consid- 
ering the limited number of lines used in this calculatibn and the 
exclusion of the H constants which could not be determined from these 
values due to its small contribution to the closely spaced Q branch 
lines. 

5. The unresolved Q branch of was observed but was too weak 
for good quantitative measurements. No S or 0 lines could be observed. 

6. Integrated intensities of the Qj (0) and S^l) lines have been 
used to calculate the polarizability constants a and y for H 2 . The 
value for y agrees with previously measured values. The value for u 
does not agree as well. The reason for the discrepancy is not known. 

7. A highly automated method for measuring the response function 
of the spectrometer has been developed using some modified techniques 
described by Jansson (16). The method uses digitized data recorded on 
magnetic tape. The Fortran coded program will average several determi- 
nations of the desired line profile and remove the calculated Doppler 
contribution by deconvolution methods. The finalized response function 
is normalized, excess base line is dropped and the results punched on 
cards for use in the deconvolution of line profiles. Calcomp plots are 
utilized for visual inspection of the progress of the program. 

8. A comparison of the response functions of two slit width settings 
are compared and show significant although small differences in profile. 
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9, A program similar in structure has been developed for using 
the response function calculated by the program above to deconvolute 
observed spectral lines. The program is again highly automated. The 
deconvoluted line is converted to an absorption coefficent curve and 
the half width and area under the curve is calculated by the Fortran 
coded program. 

10. The Q i CO) and Qj (1) , H 2 absorption coefficient line widths 
calculated by the program above have been observed as a function of 
electric field at constant pressure. The experimental error was 
larger than hoped for but definitely shows that the widths of the 
electric field induced lines observed do not depend on the electric 
field intensity. The experimental error is attributed to a poor signal 
to noise ratio. 

11. A comparison of the process of deconvolution on the Q 1 (l) H 2 
line observed at two different slit widths was made. The observed 
final widths were different by 9 percent and the integrated areas by 
7 percent. The 9 percent variation in the width is about the order 
of the experimental error encountered in the reproducibility of the 
line widths. 

12. The Q l ( 0) and Q x (l) H 2 line widths were deconvoluted and the 
widths of the absorption coefficient curves were measured as a function 
of density. Relatively few points were used in comparison to previous 
work and the extrapolated lines did not go to zero at zero pressure. 

The measured broadening coefficients did however compare favorably to 
previous results within experimental error. 
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13. A check on the behavior of the integrated area as a function 
of the electric field intensity squared at constant pressure and on the 
integrated area as a function of the product of the electric field 
intensity squared and the pressure was very satisfactory. The integrated 
areas were taken from the calculated absorption coefficient curves of 
the Q x (0) and Q x (1) , H 2 lines, and were found to vary linearly with the 
field intensity squared as well as with the product of the field 
intensity squared and the density. The intercepts at zero pressure 
and zero electric field were very close to zero area. The data obtained 
from the integrated areas were used to calculate a Lorentzian profile 
and compared to an experimentally measured profile with satisfactory 
results. 

5.2 Suggestions for Future Studies 

Two attempts to study weak electric field: induced spectra in this 
experiment (N 2 and HD) were unsuccessful. The weak N 2 ' Spectra could 
probably be seen successfully by employing the a.c. electric field 
technique used by Brannon (5) and others. In this technique a rotating 
contact is brought into the close proximity of another contact creating 
a spark gap arrangement which switches the field. If a reference signal 
is taken off. of the rotating shaft in a manner similar to that described 
for the light chopper built for this experiment then this signal can be 
used as a reference signal for the lock in amplifiers. In this way only 
the difference in the absorption with the field on and the field off 
need be amplified by the electronics and the job of offsetting the 
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background signal is eliminated. One decided disadvantage of this 
technique is that it would make intensity measurements difficult. For 
the detector employed in this laboratory, the chopping rate should 
ideally be high (1000 c.p.s.) but this would have to be balanced against 
the rise time of the chopped field. It is felt that the purchase of 
pure HD and the construction of this electric field chopper would be 
worth the value of the data obtained. 

The flushing cover used could not be completely purged of absorbing 
atmospheric gases. Since the absorptions which need to be observed in 
the HD region are small and fall in a region containing H^> absorption, 
this problem needs to be solved. The cell could be put into the 
spectrometer but the high voltage, high pressure and normal maintainance 
seem to preclude this. A more compact, source optics section in an 
evacuable chamber would solve this problem. Since N 2 lies in a region 
containing C0 2 atmospheric absorption, this would help in this area also. 
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APPENDIX A 


The program calculates several useful parameters for use in 
intensity measurements as well as frequencies for homonuclear diatomic 
molecules. The main program makes some calculations and calls on the 
various subroutines. Subroutine BOZMAN calculates the relative popula- 
tion of the various J levels discussed in Section 2.4. Subroutine 
STRONG calculates the strength factors discussed in Section 2.2. 
Subroutine FREQCY calculates the line positions in the Q, S and 0 
branches discussed in Section 2.1. The frequencies and relative 
intensities of the strongest lines are plotted by subroutine PLOTTER 
in a line printer plot. The lines are identified by their frequencies. 
The necessary input data and format are described in comment cards at 
the start of the Fortran coded program. 
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PROGRAM FOP CALCULATING LIME POSITIONS AND RELATIVE INTEN- 
SITIES. OF A HOMOMUCLEAR MOLECULE WHOSE CONSTANTS ARE KNOWN 
THE PROGRAM GIVES BOLTZMANN FACTORS, STRENGTH FACTORS, 

LINE POSITIONS AMD RELATIVE INTENSITIES FOR. THE FIRST 26 
LINES IN EACH BRANCH FOR BOTH PARALLEL ANC PERPENDICULAR 
POLARIZED RADIATION. THE STRONGEST LINES ARE PLOTTED OUT 
AS FREQUENCY VS RELATIVE INTENSITY. ORDER AND DETAIL OF 


DATA CARDS IS 
CARO 1 

CARD 2 
CARD 3 


AS FOLLOWS. 

...TITLE CARD FOR THE RUN, TYPED ANYWHERE 
ON CARD 

. • . BV IN COLUMNS 1 TO 15 MOLECULAR CON- 
STANT 3 F°R GROUND STATE V 
DV IN COLUMNS 15 TO 30 MOLECULAR CON- 
STANT D FOR GROUND STATE V 
HV IN COLUMNS 30 TO 45 MOLECULAR CON- 
STANT H FOR GROUND STATE V 
...T IN COLUMNS 1 TO 15 KELVIN TEMPERA- 
TURE 

SO IN COLUMNS 15 TO 30 SP IN DEGENERACY 


FOR COO J LEVELS 

SE IN COLUMNS 30 TO 45 SPIN DEGENERACY 
FOR EVEN J LEVELS 

CARD 4....AV IN COLUMNS 1 TO 15 AVERAGE POLARI- 


CARD 5 

) 


CARD 6 



Z1BILITY 





GV 

IN COLUMNS 

1 

TO 

15 

ANISOTROPY OF 


POLARIZATION 




VO I 

IN COLUMNS 
STANT 

1 

TO 

15 

VIBRATIONAL C ON- 

BVP 

IN COLUMNS 

15 

TO 

30 

MOLECULAP CON- 


STANT B CF 

EXCITED 

STATE V PRIME 

BV 

IN COLUMNS 

30 

TO 

45 

MOLECULAR CON- 


STANT B OF 

OR' 

3 UNO 

STATE V 

DVP 

IN COLUMNS 

45 

TO 

60 

MOLECULAR CON- 


STANT D OF EXCITED STATE V PRIME 
DV IN COLUMNS 1 TO 15 MOLECULAR CON- 
STANT D C c GROUND STATE V 
HVP IN COLUMNS 15 TO 30 MOLECULAR CON- 


STANT H OF EXCITED STATE V PRIME 
HV IN COLUMNS 30 TO 45 MOLECULAR CON- 
STANT H CF GROUND STATE V 
CARD 7.... SC ALE OF LINE PLOT IN 15. SUBTRACT THE 
ROUNDED OFF FREQUENCIES OF THE TWO 
CLOSEST LINES. THE SCALE FACTOR TIMES 
THIS NUMBER SHOULD BE AT LEAST EQUAL 
TO ONE 


DCUOL EPREC I SI ON F ( ?6 ) , S 1 10 { 26 ) , SLQ ( 26 ) , S 11 S ( 26) , 
1SLS (26 J ,$ 110 (26) , SLCM 26) , V0( 26) , VS { 26) , VO l 26) , 

1A 11Q( 26 ) , ALQl 26) , A1 IS (26 > , ALS t?6 ) , A1 10(26) , ALD126 ) , 
IN A 1 IQ { 2 6 ) » fJ ALQ ( 26 ) , N A 1 1 S ( 2 6 ) ,N AL S l 26 ) , N A 1 1.0 ( 26 ) , 

IN ALC (26) ,MAX(73)*AA(73),VV(78)»VVV(78)»AAA(73)» 

1BIG (6) ,HUG E 



INTEGER IBAX, IFLAG, I AV(?3) 

CALL BCZMANIF) 

CALL STRONG (SI 1C,SLC,S11S,SLS, SI lOtSLO) 

CALL FFcGCY(VO*\JS t VC) 

Vi PI TE C 6 » 1 CO J 

100 F C R ’ 1 A T ( 1 H 1 ♦ 23HBCLTZ V ANN*STRENGTH*FPEQUENCY> 

vs rite it, rci ) 

101 FCRVAT (1HC,20X,3HQ BRANCH » 2 OX ,8HS BRANCH, 20X, 

1 8H0 BRANCH) 

fcRITE (6 ,102 ) 

102 FORMAT ( 1H0 ,7X 8HPARALLEL,4X, 13 HP ERPENOICULAR 
14X,8HP A RAIL EL ,4X, 1 3 HP ERP EN C ICIJL AR , 4X , 8HPARALLE L , 4X , 
113HPEFPENCICULAR) 

CC2J = 1 » 26 
I =J-1 

A110tJ.) = F{ J)*S110C J)*VQ( J) 

ALC(J)=f (J}*SL0<J)*VQ(J) 

AHS(J)=F<J)*$US(J)*VS(J) 

ALS{J)= C {J}*SLS{J)*VS(J) 

AUCtJ) = F{ J)*S11C( J)*VC( J) 

ALO(J)*F(J)*SLG{Ji*VQ(J> 

2 WRITE (6, 106) I , A11GI J) , A LG ( J > , All S t J ) , AL S ( J ) , A 1 10 ( J ) 
1 ALC ( J ) 

106 FCRMAT { 1B0 , 15 ,6015 .6 ) 

N = 1 
N =26 

CALL .VAXOH AUQ,«, IMAX} 

BIG(N)sAllOUMAX) 

N = N + 1 

CALL i-’AXOK ALQ, v , I MAX) 

6 IG(N) = ALC( IM AX ) 

N=N+ 1 

CALL MAX01(A11S,M, IMAX) 

BIG IN) =A1 IS ( I MAX) 

N = N + 1 

CALL f-'AXOl ( ALS ,M, I^AX ) 

B ! G { N ) = A L S ( I ,J! A X ) 

N = N + 1 

CALL VAX01 ( A110, M , INAX) 

B I G t N > — A 1 1Q( I M A X ) 

N = N+ 1 

CALL MAXO 11 ALH, M, I MAX) 

B I G ( N ) = A L 0 ( IM AX ) 

M = 6 

CALL MAX01(BIG,M, I MAX) 

H L C- E = 5 1 G ( I OAX) 

C C 5 J= 1 , 26 

NAUQU) = A11Q(J)/HUC5 
K ALC { J) =ALO (J) /HUGE 
NA11SIJ}=A11S{J) /HUGE 
NALSt J)=ALS(J )/HUGE 
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NA110! J)=A110! JJ/HUGE 
5 NALGC J) =ALOl Ji/HUGE 
WRITE<6,200) 

200 FORMAT C 1H1» 31 HNORMAL I ZED PROOUCT OF B,S AND Vi 
WRITE (6, 1 Cl) 

WRITE! 6 » 102 ) 

DC 11 J = 1 » 26 
I = J-1 

11 WRITE 16, 106 II ,NA11GIJI»NALQ! J I, NAilS! J) »NALS! J), 
1NA110! J) tNALOt J) 

CALL OROER<NAlIQ»NAllS.NAilO*VQ*VS, VG,VV,AA) 

J^l 

K=78 

301 IF( J.EQ.79) GOT 03 00 
CALL MAXO 1{ WtMflMAX) 

IF ! VV MM AX) •EQ.0.0D0JGOTO300 
IF! IMAX.GT. 26IG0TQ600 
GOT 0601 

600 IF!IMAX.GT*52) GOT0602 
GGTC6Q3 

602 IF LAG =3 
GGT0604 

601 IFL AG=1 
GCTC604 

603 IFLAG=2 

604 I J) “I FLAG 
VVV!J)=VV! IMAX) 

VV( IMAX)=0.0D0 
AAAtJ) =AAU«AX) 

J=J + 1 
GCT0301 

300 K= J - i 

WRITE (6 ,304) 

304 FORMAT 1 1H1, 19HLINES TO BE PLOTTED) 

WRI T E 16 » 305 ) 

305 FORMAT! 1HG»2X»4HL I NE,6X, 9HFREQUENCY » 16X »9HI NTENS IT Y ) 
DO 302 J=1»M 

302 WRITE(6,303) J,VW! J) ,AAA{J) 

303 FORMAT! IHQ, 15,2020.9) 

WRITE 16,1100) 

1100 FORMAT! IHii 31HRELATIVE INTENSITY VS FREQUENCY) 

WRITE (6,1101) 

1101 FORMAT! 1HC,38H FOR PARALLEL POLARIZED RADIATION) 

CALL PLOT £R ! AAA* VVV » I AM »M ) 

STOP 

END 
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SUBROUTINE 80 Z M A M { F ) 

DCU8LE PRECISION FV J , BOLT ,NUM HOD t SUM, SP IN, NAME! 10 ) ,X» 
1CCAL.FI26) 

P E AC (5,E0C)NAME 

800 FCRVAT(IOAC) 

WRITE! 6, 801 INANE 

801 FORMAT (1H1 »35X, 10A8) 

ViPITE{fc,30?) 

3 02 F CFMAT ( 1H0 » 35X , 3 5 H*** ****$********* ******** ********** ) 
7 CC PEAC(5i50C)BV,nVtHV 
500 FORMAT (3 E 15.0) 

WPI T E < 6 » 3 C 0 ) 

300 FORMAT (1HC» 15HINPUT DATA t BV=, 20X.3H0 V= » 20X ,3HHV= ) 
WRITE(6,301 )3V,DV,HV 

301 FORMAT { IHCt 3025.9) 

RE AC { 5 » 500 ) T t S 0 1 S E 
U P I TE ( 6 , 600 ) 

6C0 FORMAT! 1H0, 14HINPUT DATA, T= , 20X , 3H SC= » 20X , 3HS£ = ) 

WRI TE (6 ,301)T,SC,SE 
WPITE{fc,40C) 

400 FORMAT ( 1H0 , 4X , l 8J , 14X , 1 7HB0LT ZMANN FACTORS) 

J = 0 
IC0=1 
5 K= J+ 1 

F V J = J*K *8V- J**2*K** 2* DV + J* * 3*K**3*H V 

CC AL= ( 1 . 3 80 622 0-16 * T ) / { 6 . 62 1960-2 7*2 .99 7925D+ 10) 

X=FVJ/CCAL 

BCLT=DEXP(-X) 

13 DEG=2*J+1 

GOTO (10,15) , I GO 

10 I GC = 2 

S P I N = S E 

11 J=J+1 

IF( J.E0.102 JGCT020 
NUM( J)=e 0 LT*DBG*$P1M 
I F ( N U M ( J ) .EC.O.ODO )GCTC5 
IF (NUfM J) .LT. 1. CD-20) GOTO? 

C-CTC5 

2 MUM ( J ) = C . GOO 

J=J + 1 

I F ( J.E0.102 ) G G T C 2 0 
GCTC2 

15 S P IN = SC 
IGC = 1 
CCTC11 
20 SU^=0 

CC25J-I , 101 
25 SUR=SUV+MUM( J) 

0C30J= 1 , 26 

F (J )=NUMl J ) /SUM 

I=J-l 



30 W R I T E < 6 » 5 03 } I » F C J 1 
503 FCRMAT(1H0,I5»D25.6) 
RETURN 
END 
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SUBROUTINE STPOMG(SllC»SLC»SllStSLSt SllO.SLC) 
OOUBLEPRECt SION S 1 10 ( 26) , S LQ t 26 ) , S 1 1 S ( 26 ) , SI $ ( 26 ) , 

IS 110(26) , $ L G ( 26 ) ,A,9,C,C,E,F f G,H,Q,GV,R»S,T ,U,V,W»AV, 
lPtO 

READ (5, 3 CO) AV,CV 

300 FCRPATI2E15.0) 

WRITE (6 ,301 ) 

301 FORMAT { l H 1 » L 5 H T U P UT DATA, AV= ,20X,3 HGV= ) 
WRITE(6,302)AV,GV 

302 FCRMAT ( 1H0, 2025.9) 

WRITE (6, 303) 

303 FGRVAT ( 1H0 , 35X , 16HSTRENGTH FACTORS) 

J = 0 

R=Q .000 

18 S=R+ 1 . OCO 
7=2 *R- 1 .000 
U =2* R+ 3 . 0 DO 
V=4.0D0/45. GDO 
W=1 .0C0/15.0D0 
J = J + 1 

S11QI J )=l R*S) /(T*U)*V*(GV#*2)+IAV**2) 

SLQI J) = ( R * S ) / <T*U)*Vi*tGV**2 ) 

IF{ J.E0.26IG0T01S 
R=R+1. OCO 
G C TO 18 

19 J = 0 
A-O.GDO 

P=2.0Q0/1 5.0D0 
Q = 0. 1CDC 
5 B=A+1«0C9 
C = A + 2 . 0 GO 
D=2.0DC*A+3.0C0 
E=2.0C0*A+1 .OCO 
J = J + 1 

SI IS { J ) =( 3*C )/ ( D*E »*P* tCV**2 ) 

S l 5 (J > = ( 3*C > / { D * E ) * C + ( G V# *2 ) 

IF{ J.EC.26IGGT012 

A-A+1 .OCO 

GCTC5 

12 S11C (l ) =0.0 CO 
stem =C.ODO 
J = 1 

F = 1 . OCO 
15 G=F-1.CC0 
H=2*F-1 .000 
r=2*F+ 1.000 
J = J + 1 

S110 IJ) = 1 F + G ) / ( H* G ) * P * ( G V ^ * 2 ) 

SLG ( J } = { F *G ) / ( H *0 ) * Q * ( G V * * 2 ) 

1 F { J.EC.26) G0T016 
F=F+1.0C0 



GOTO 15 

16 WRI TE 1 6 » 504) 

504 FORMAT i IHO* 20X* 8HQ BRANCH* 20X* 8HS BRANCH*20X» 

18HC BRANCH I 

WRI T£ ( 6* 505 ) 

505 FORMAT ( 1H0,3X,2HJ,4X,8HPARALLEL *6X, 1 3HP ER PENDICULA R 
1 4X» 8HPARALLEL , 6X , 13HPERPENDICULAR * 4X * 8HP ARAILEL* 6X » 
113HPERPENDICULAR) 

DO 17 J = l * 26 
I = J»1 

17 WRITE (6*5 06 il»SllQiJ)*$LQ(J)*SllS4J)*SLSCJ)»S110tJ) 

1SL0CJI 

506 FORMAT ( 1HQ* I 5 *6015.6) 

RETURN 

END 
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SUBROUTINE FREOCY(VC»VSfVO) 

DCUSIEPRECISICN VO f 26 ) , VS ( 2 6 ) , VO ( 26 ) , V < 26 ) , A, B, C , D 
REAC (5, 700) VO I , RVP , BV,DVP 
7 CO FCRVAT (4E15.0) 

• PEACI5, 701) DV,HVP,HV 
701 F C F V A T ( 3 E 1 5 • 0 } 

WRITE ( 6 1 3 C 0 ) 

300 FCRtYA T { 1H l, 16 HI NPiJT DATA, VD I = , 20X , 4HBVP= , 20X , 3HB V=) 
‘WRITE (6 ,301 ) VO I ,BVP,BV 

30 1 FORMAT ( 1H0, 3D25.5) 

WRITE (6,302) 

302 FORMAT ( 1H0, 16HINPUT DATA, Q VP = , ?QX , 3HDV= , 20X , 4HHVP= , 
1 20X , 3FF V- ) 

kRI T E { 6 , 3 03 ) :) VP , D V , K V P » H V 

303 FORMAT ( 1H0, 4025.9) 

WPI TE (6 ,3 34) 

304 FCF^AT ( IHGf 35X, 14HL INE FREQUENCY) 

J = 0 

A = C*CDC 
J GC = 1 
I DC\E= 1 
5 6 = A 

C*A«-1.0C0 
D=A+ l, 000 
11 J=J+1 

IF( J.EQ.27)G0TC16 

V ( J ) =V 0 1 + 3*0*'3VP-A*C*eV-5*< ! 2*D**2*DVP + A*< ! 2*C**2*DV+ 
16** 3 * D ^ 5 - 3 * H V P - A * * 3 * C * * 3 * H V 

A = A + 1 . C C 0 
G CT C ( 5 , 13 ,25) , I GC 

16 GOTO ( 17,15,30) , IPONE 

17 C C4Q J = 1 ,26 

40 VC(J)-V(J) 

J = 0 

A ♦ 0 CO 
IGC = 2 
I CC KE = 2 

18 C=A+1.CC0 
R=A+2.CC0 
C - A 4- 3 . 0 C 0 
GCTCll 

19 CE41J=1,26 

41 VS( J)=V(J) 

V { 1 JsO.CDO 

V ( 2 ) =0. ODD 
J = 2 

A =2 . OCO 
I GC = 3 
I CC \E=3 
25 C-A-l.CDO 
e-A-2,0C0 



D-A-3.0D0 
GOTO 11 

30 D042J= l *26 
42 V0IJ)=ViJ) 

WRITE (6,507) 

507 FORMAT! 1H0 , 1HJ,7X»8HQ BRANCH, 12X»8HS BRANCH, 12X 
18HG 8RANCH ) 

CC3 1J= 1 ,26 
I=J-1 

31 WRITE (6,503JI ,VQU) »VSC J) ,V0( J) 

503 FORMAT ( 1H0, 1 5 ,3D2 0»9) 

RETURN 

END 



SUBROUTINE MAXOl CA,M, INAXI 

DOUBLE PRECISION GREAT, A{73> 

INTEGER IMAX 

GREAT=C.ODO 

I MAX=1 

DO 1 J=1,M 

IF (AC J)*EQ.O«ODOJGQTOL 
I F I A< J ) *LT. GREAT) GGTG1 
IMA X= J 
GREAT»A ( J ) 

1 CONTINUE 
RETURN 
END 
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SUB ROUT I NE QRDE R INA1 1C, NAl 1 S » NA1 10, VQ , VS » VO, VV f A A ) 
DOUBLE PRECISION NA11QI 26) *NA11S( 26) ,NA llOi 26} ,VQI 26) » 
1VSC26) 9 VO ( 2 6 ) ,VVC781 ,AAC78) 

J=1 

DC 3 1=1,26 

vvtn=vei Ji 

AAt I)=NA11QIJ ) 

3 J-J+l 
J = 1 

DC 4 1=27,52 
VVU)=VS(J) 

AA( Il-NAllS(J) 

4 J=J + 1 
J=1 

DC 5 1=53,78 
VV(I)=VOI Ji 
A A ( I l=NAI10< J ) 

5 J=J + 1 
1 = 1 

8 IF! AAC I). EQ.0.0D0IGGT06 
IF! A AC I).LT„0,10-3)G0TQ6 
GOT 07 

6 VVI I )=0»0DQ 
AA( I 1=0.000 

7 1=1+1 

!F< I.EQ.79IG0T09 
6 CTO 8 

9 RETURN 
END 



SUBROUTINE PLCTER ( A’, V , I AM, M> 

DOUBLE PRECIS ION V ( 76 ) , A { 78 ) , FREQ , R, P'O, SC ALE , F I RST 
INTEGER I FREQ, IPQ, J, I , ISTEP, LINE * 101 J , IAM{?8) 

READ ( 5,400) SCALE 
400 FORMAT (1EI5.0) 

' WRI TE (6,406) 

406 FORMAT ( 1H0, 6HSC ALE= 3 
WRI TE (6 ,402 JSCALE 
402 FORMAT { 1KO,D20. 9) 

CAT A IPL, IAST, TMX/2H ,2H* , 2HX / 

DC91CJ =1,101 
9 10 LINE (J 3 = 1 AST 
921 ViPI TEC 6, 912) LINE 

912 FORMAT ( 1H0, 101A1) 

I S7 E P = 1 

M*M+l 
J = 1 
1 = 1 

801 IF (V ( J ) .EQ.O.0DO3GCTDSCC 
GCTC91 1 
800 J=J + 1 
1 = 1 + 1 
GCT0801 

911 FIPST=V(J) 

122 R=DARS(V(n -FIRST) 

F R 5 Q = ( R * S C A L E 3+0.55D+C1 
I FREQ=FPEC 

PC= A ( I 3*0. 100+03+0. 1500+01 
I F G = P G 

DC9 19 J=l,101 
919 L IN E ( J ) = I OL 
L I NE { 11=1 AST 

300 IF{ IFREG.EO.l STEP JGCT0913 
GCTC913 

913 DOS 14 J = 1 , 1 P Q 

914 LINE(J)=l M X 

IF (I AM ( I) .50. 13 GOTO 30 1 
I F ( I AM ( I ) .50.2 3G0TQ303 
IF (I AMU ) .EG. 3 3 GOT 0305 

301 WR I T 5 ( 6 » 3C23L INE,V( I ) 

302 FORMAT UH*, 101AI ,2X , IRQ, D15.6) 

GO TO 30 7 

303 VsRITE(6,304)LINE,V( I) 

304 FORMAT (IH , 10 1A 1 , 2X , 1 HS ,01 5 .6 3 
GOTO 30 7 

3 05 IvR I TE ( 6 , 3 0 6 3 L I NE , V ( I 3 

306 FORMAT ( IH , 101 A 1 , 2 X , 1H0 ,D15.6 ) 

307 1=1+1 
IFU.EO.M3GOTP917 
ISTEP=ISTEP+1 

I F I ISTEP.E3.400 3GOTC917 



G0T0122 

918 WRI TE(6t9l6jLINE 

916 FORMAT { 1H , 101AU 
I STEP=I STEP+1 
IFCISTEP. E 0.400 I G0TG917 
GOT 0300 

917 RETURN 
END 



APPENDIX B 


The program for the response function is coded in Fortran and 
runs on the IBM 360/65. Two minutes are required with 13 data runs 
and a total of 15 iterations. Three passes in loop 1, five passes in 
loop 2, and seven passes in loop 3 are used. The main program manip- 
ulates the data and calls on the various subroutines. The program 
first calls the subroutine READER which reads the first and the second 
specified base line block for the first data set. These points are 
averaged and divided by 1024 since the digitized data is recorded as 
values from 0 to 1024. This determines the base line for the first 
data set which is then read in by the subroutine READER. The base line 
value is subtracted from the data and the input value of the 100 percent 
absorbing point is used to convert the data to a fractional absorption. 
The data is smoothed once using the subroutine SMTHER which is part of 
Janssen's (19) deconvolution program. The array is then searched for 
the maximum value which is the peak absorption and the data points 
corresponding to the half absorption values are then found. This data 
is used to calculate the line center. The next steps in the program 
align the several data runs, with respect to this first computed line 
center point. The first data set is not moved but subsequent data sets 
are aligned so that the measured center point of the line lies at this 
first computed center point. Two separate loops, one for moving the 
data to the left and one for moving the data to the right, are used. The 
results for each subsequent data set is added into the array D, for 


181 



182 


averaging after all data sets are processed. The subroutine BLOCK is 
called three times in these steps and plots the raw data, the smoothed 
data and the aligned data for each data set. After looping back to 
read each subsequent data set the array D is averaged to give an average 
line profile. Subroutine GAUSSIAN is used to calculate the doppler 
profile for use in the deconvolution steps. The doppler profile is 
normalized and used in calling the subroutines DECONER and DECONVO which 
are part of Jamison’ s program which have been modified as discussed 
in Section 4.7. The half-width of the deconvoluted profile is measured 
by looping back into the first part of the program which finds the points 
corresponding to the 50 percent peak absorbing points. All excess base 
line is then discarded which lies outside four times this computed half 
width. The area under the curve is then computed and each point normal- 
ized to this value so that the area is equal to unity. The final result 
which is the response function is punched on cards for use in the 
deconvolution program. Subroutine BBPLOT is called twice in these 
final steps. This subroutines plots the averaged data and the deconvo- 
luted data with the excess base line dropped. It may be necessary to 
drop a few more base line points before using the response function in 
the deconvolution program. The data needed for the program and the 
format of this data is documented in comment cards listed at the start 


of the program. 
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C THE PROGRAM TAKES UP TO 13 RESPONSE FUNCTIONS FROM 
C MAGNETIC TAPE AMD AVERAGES THEM AND CECCNVCLUTES THE 
C DOPPLER WIDTH FROM THE AVERAGE. THE PROGRAM FIRST 
C READS ID, THE 100* ABSORBING POINT. FOR THE FIRST 
C RESPONSE FUNCTION, TWO S°ECIF1ED BLOCKS ARE READ AND 
C AVERAGED. THE INPUT DATA IS READ AND SMOOTHED ONCE. 

C THE BASE LINE IS SUBTRACTED AND THE DATA IS CONVERTED 
C TO A FRACTIONAL ABSORBANCE. SUCESSIVE INPUT DATA IS 
C READ AND THE RESPONSE FUNCTIONS ARE ALIGNED, AVERAGED, 

C CECCNVCLUTEO, NORMALIZED AND PUNCHED ON CARDS. THE 
C WIDTH OF THE RESPONSE FUNCTION IS CUT TC 4 HALF 
C WICTFS IN TFE FINAL STEPS. CALCCMP PLOTS ARE MADE AT 
C INTERMEDIATE STEPS TO CHECK THE PROGRESS OF THE 
C PROGRAM. 

C DATA CARDS AS FOLLOWS: 

C CARD 1, LABEL FOR RUN AN V WHERE CN CARD 
C CARD 2, NUMBER OF RESPONSE FUNCTIONS' TO AVERAGE 
C CARD 3, NUMBER OF BLOCKS READ FOR I G I = 1 ) IN 5, 

C STARTING BLOCK IN 10 

C CARD 4, NUMBER OF BLOCKS READ FOR FIRST BASE BLOCK 
C ( = 1) IN 5, STARTING BL.CCK IN 10 

C CARC 5, SAME AS 3 FOR SECOND BASE BLOCK 

C CARD 6, NUMBER OF BLOCKS READ FOR RESPONSE FUNCTIONS) 
C IN 5, STARTING BLOCK IN 10 

C NEXT CARD, REPEAT 4 TC 6 FOR EACH RESPONSE FUNCTION 
C NEXT CARD, tf OF CHARACTERS IN FIRST BB P LOT LABEL IN 5 

C NEXT CARD, LABEL FCR ABOVE STARTING IN ZERO 

C NEXT C A p D » DATA TO CALCULATE GAUSSIAN 
C FREC. OF LINE IN 15, *-■ ASS OF MOLECULE IN 30, 

C KELVIN TEMP. IN 45, FREC. PER DP. IN 60 

C NEXT CARD, DECONVOLUTION DATA 

C NEXT CARD, # OF CHARACTERS IN SECOND B8PL0T LABEL IN 5 

C NEXT CARC, LABEL FOR ABOVE STARTING IN ZERO 

C NEXT CAPO, # OF CHARACTERS IN THIRD 8RPL0T LABEL IN 5 

C NEXT CARD, LABEL FOR ABOVE STARTING IN ZERO 

R E A L * 4 10 

COMMON /A/X U 280 ) ,J,SRMS/B/Y{ 1280) ,A( 150) ,NA 
DIMENSION IBUFI 664) , 0{ 1280 , NAME (20) 

INTEGER I MIDI 15) 

CALL PLOTS! I BUF ,664,80.0, ’BOYDS * , 'PEMGTE2#' ) 

CALL PLOT (0.0, 0.5, -3) 

PE AC ( 5, 7 23) NAME 

723 FORMAT (20A4) 

WPI TE (6, 724) NAME 

724 FORMAT ( 1H1, 20A4) 

900 F C F M A T ( 1H , 10F12 .5 ) 

132 FCPMA1 (215) 

618 FORMAT ( 1 HO , 15X , 3HNB = , 13, 10X,4HNB1=, 13) 

725 FORMAT { 1H0 , 2 5H - Tr **** * * # a « :::****;* ** ) 

DC 800 1= 1, 1280 

300 C ( I ) = 0 .0 
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PEAC (5, 139) ITEST 
139 FORMAT ( 15) 

C READ ic 

R E A C ( 5 » .1 3 2 ) N B ♦ N B 1 

WP1 TF. (6t613)NB,N31 

CALL READER l N B » N3 1 ♦ G » 815) 

TCTAL=Q.O 
DC 702 I = 1,256 
702 TCT AL = TCT AL+X ( T ) 

IC=TGTAl/<0.256E+03*0.1024E+04) 

WRITE (6, 5 02 ) I C: 

V, R I T E ( fc » 7 2 5 ) 

502 FCPMAT(lH0t2OXt3HIf>,FlC.6) 

C READ BASE LINE FROM TWO SEPARATE BLOCKS 
I A M s 1 
714 ICC=1 

B ASE=Q. C 

245 PEAU5, 132)N8»N8l 
WRI TE (6 ,618) NBt N n l 
CALL READER (NB»N8l»C»£15) 

DC 712 I -• 1 » 2 5 6 
712 BASE=BASP+XtI ) 

I F ( I GO -EG .2 ) GO TO 247 
I GG = 2 
CC TC 245 

247 6ASE=BA$E/{ 512.0*1024.0) 

WRITE! 6 t 833)3ASE 

83 8 F CF -7 AT ( 1H 0 » 20X T 5HB A SE = , FI 0.6) 

C READ INPUT DATA 

REAL ( 5 » 1 3 2 ) N 3 * N 5 1 
J=NB*256 

WRITE! 6,6 1 9 ) N B » N 3 1 

619 FORMAT! 1H0, I 5X , 3 H N B - » 1 3 , 1 OX f 4HNB 1 = t 13 ) 

ViP. ITE (6,7C1 ) 

7 Oi FORMAT ( 1H0 , I7HI NPUT INTENSITIES) 

CALL PeAQEP{*!fcitN31fC,Sl5) 

DC 938 1= I * J 

933 XII) =X( I) /0.1C24F+C4 
C PLCT RAW CAT A 

CALL PBLCCK(XtJ) 

CALL PLCT16. 0,0. 0,-3) 

C SMOOTH INPLT DATA 
CALL S M T H E R ! 9 5 ) 

C SUBTRACT EASE LIME AMC CONVERT TO FRACTIONAL ABSORPTION 
R I C = I 0 - E A S E 
370 CC 788 1= 1, J 
788 X! I ) = ! X! I ) -EASE ) /P IC 
C WRITE INPUT INTENSITIES 

901 WRITE (5 ,9 00 ) IX ! I ) » 1 = 1 » J, 10 > 

C PLCT SMOOTHED DATA 
CALL B8L n CK!X,J) 



185 


CALL PLCT (6 .OtO.Ot-3) 

C ' F INC MID POINT 
3 £ 5 GPEAT = C.O 

CC 713 1=1, J 

I F ( X ( I ) .LT.GPEATIGO TC 718 
GPEAT=X (I ) 

ICRGAT= I 
713 CENTIME 

IF ( I GO .53 .4 )Gj TO £64 
HALF =G RE AT/2.0 
I*IGPEAT 

768 1=1+1 

IF(XU)-HALF) 770,770,763 

770 IHALF2=I 
I =1 GREAT 

769 1=1-1 

I F t X ( I ) -H AL F ) 771,771,7 69 

771 rhALFl=I 

IKICUAM)=<I HAL F1 + IFALF2 ) / 2 

WRITE ( 6,772 ) I GR EAT, I HAL FI , IHALF2 , I MI Ot I AM) 

772 FOR MAT ( 1H0 , 7HI GR E AT= , I 9 , 10X , 7H I HALF 1 = , I 5 , 10X, 

1 7HIHALF2 = , 17, 1 CX , 5H I M I C= 1 7 ) 

HRI TE (6,7 25) 

IR£00 = I V IC(1)-I V 1D{ I A >V ) 

888 IF t I RE DC > 70,71,72 

70 CC 969 1=1, J 

k=i-irecc- 

IFIK.LT .1. OR. K.GT.J )G0 TO 697 

xu )=xm 

GC TO 3 69 
597 X(I )=C.OEQ 

969 CCNTINLE 
GC TC 71 

72 DO 693 L= 1 , J 
I = J -L + 1 
K= I-IRRCC 

IFIK.LT.l.OR.K.GT. J)GO TO 699 
XU )=X (K) 

GC TC 698 
699 X ( I } = 0 . OEO 
698 CONTINUE 

7 1 I F ( I GO . E9 . 3 ) GO TO 999 
C PLCT ALIGNED DATA 

CALL 6 B LOCK ( X , J ) 

CALL PLOT (-12. 0,0. 5, -3) 

C ADD RESPONSE FUNCTIONS 
CC 9 70 1= 1, J 

970 c (i )=cm+xm 

1 F II AM . SQ. I TE ST )G0 TC 713 
I A N = I A R + 1 

C RESTART FOR NEXT RESPONSE FUNCTION 
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C-0 TO 714 

C AVERAGE RESPONSE FUNCTION 

713 CALL PLOT (18.0, -FLOAT *ITE$ 1—21*0. 5,-3) 
CUTEST 
DC 971 I=i, J 
971 XUI=OIII/C 
WRITE 16,782) 

782 FORMAT* 1HG, 26HAVER AGED RESPONSE FUNCTION) 
WRITE (6 ,900) CXt I), 1*1, J, 101 
C PLOT RESPONSE FUNCTION 
CALL 88PLCT (256.0) 

C CALCULATE GAUSSIAN FOR DECONVOLUTION 
CALL GAUSS* A, NA) 

C WRITE GUT GAUSSIAN 
WRITE *6, 7261 

726 FORMAT * 1H0* 15HD0PPL ER PROFILE) 

WRITE * 6, 9001 1 A* I)»I=1»NA) 

C OECONVOLUTE RESPONSE FUNCTION 
CALL DC ON VO 

C PLOT DECONVOLUTED RESPONSE FUNCTION 
CALL BB PLOT (256.0) 

C DROP EXCESS BASE LINE 

L ENGTH=2* * I HALF2-I HAL FI ) 

IREDQ=-* I MID* I) -LENGTH) 

I G0=3 
GO TO 888 

C NORMALIZE RESPONSE FUNCTION 
999 ARE A=0. 0 

J=2*LENGTH 
DO 778 1-1 , J 
7 78 AREA=AR£A+X 1 1 ) 

DC 779 1=1, J 

779 X*I )=XC I) /AREA 
WRITE* 6,780) 

780 FORMAT * 1H1»23HF INAL RESPONSE FUNCTION) 
WRITE *6, 7811* XU), 1*1* J I 

781 FORMAT * 1H ,10F12.5J 
C PUNCH RESPONSE FUNCTION 

WRITE *7, 777 MX* I),I=i,J) 

777 FORMAT *8E 10.4) 

C CHANGE SCALE AND REPLOT 
I GC=4 
GC TO 885 
884 DC 886 1=1, J 
886 X* I )=Xt I) /GREAT 

CALL B8PLGT* 256.0) 

CALL PLCT*0. 0,0. 0,999} 

15 CALL EXIT 
END 
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S U6RCUT IN E R E A D £R ( N X , N8 1 , K00 E , * ) 

KX BLOCKS OF DATA FROM UNIT 9 (TAPE) — 256 1*2 FORM 

RgAO INTO X STARTING WITH BLOCK # NBl 

1 F KOBE NE 0 UNIT 9 ENDED 

-tOC = NEXT BLOCK ON TAPE, #0 = LABEL 

COMMON /A/XI 12801, KKK.SRMS 
I NT EGER* 2 INPUT! 256 I ,U 
LOGICAL*! LABEL! 141, LK,L! 21 
EQUIVALENCE !U1),LI1 
DATA LOC/C/ 

IFiLCC.EQ.O) GO TO 6 
IFtNBl-LGCJ 3,2,1 

3 REWIND 9 
L OC ~ 0 

6 IC=1 

R E AC !9 » 101 * END=7 ,ERR=71 LABEL 
L CC = l 

100 FORMAT 1128 1 2A2) 1 

101 FORMAT (14A1) 

1 NC=NB 1-LOC 
ID =2 

DO 5 1*1, ND 

LCC=LCC+1 

5 READI9, 100, END=7 ,ERR=7) M 

2 CONTINUE 
NXX=~256 
I C= 2 

CO 4 1=1, NX 

NXX-NXX+256 

READ !9 , 100, END= 7, ERR = 7 ) INPUT 
LCC=L0C+1 
CC 4 J •=l,.256 
U = INPUTU) 

LK=Ll 1) 

LI11=LC2I 

L(2J=LK 

4 X (NXX+J )=L I 

I FI KODE . E Q» 0 1 RETURN 

C -END FILE 

I D=4 

8 REA0I9, 100, END=9 ,ERR=7I H 
LCC=LCC+i 

GO TO 8 

9 LCC=0 
RETURN 

7 WRITE (6,200) ID,LCC 

200 FORMAT I5X3HID= 1 3, 5X4HL0C= 15 , 5X, 

1 ‘READER ERROR PROGRAM REROUTED*) 

RETURN 1 
END 



SUBROUTINE GAUSS IA»NA) 

DI MENSION A 1 1 50) 

READ(5,331)VO,W»T f FPERDP 

331 FORMAT (4E 15.0) 

€=2.9979256+10 
8=1.386226-16 
WRITE! 6,330) 

330 FORMAT (1H1, 23HINPUT DATA FOR GAUSSIAN) 

WRITE i 6 ,332 ) VO » W ,T , FPERDP 

332 FORMAT 1 1H0,6HFREQ.= ,E 15.7, 5X,5HMASS= »E15.7,5X 

1 6HTEMP . = »E15. 7 »5X » 1 QHFREQ ./DP. = , E15.7I 

Q={ W*C*C)/I 2.*V0*V0*B*T) #FPERDP*FPERDP 
R =- A L OG 1 0 . 1 E-3 ) 

IMID=SQRT !R) /SORT! Q) + 1.0 

NA=2*IMID-1 

S=S0RT<Q/3,1415) 

A C IMIDI =S 
MI D=I MI C— 1 
DO 1 1=1, MID 
J = I M 1 0+ I 
K=I MID-I 
X= I 

V=EXPf-X* X* 0 ) * S 
A(JI=Y 
1 AIK)=Y 
RETURN 
END 



SUBROUTINE SMTHER I NSM ) 

COMMON /A/DATAI 12801 , NO t SRM $ 
DIMENSION S (99 ) » T{ 200 1 
DATA KSW/O/ 

7 IF ( NSM. EQ «K SM ) GOTO 8 
K-(NSM-il/2 

R=NSM 
K3=NSM+ 1 

A2=K*(K+i)*<2*K+l) /6 
A=3*K# C K+ II - 1 
XNGR«=R*A-10.*A2 
DC I 1=1, K 
P = I 

1 $tli={ A-5.*R*R> /XNGRM 
SO=A/XNCRM 

K1=K+1 

8 SK=0. 

CO 2 1=1, NS M 

2 T ( I >=DATA{ I ) 

I =K 

3 1=1+1 

T<NSM+li=OATAU+Kl) 

ST=SO*TtKl I 
DO 5 1*1, K 

5 ST=ST+S (L)*CT<K 1+0 +T { Kl-L ) ) 
DO 6 L= 1 , NSM 

6 T ILI=T < U + 1 ) 

S K= S K+ f 0 A T A 1 1 J ~ S T J * * 2 
OATAM) = ST 

IFC I+Kl .LT.ND) GO TO 3 
Z=I-K 

SRMS=SQRT C SK/Z I 

KSH=NSH 

RETURN 

END 



SUBROUTINE BBLQOK I X , N X) 

DIMENS ICN XC1)»S{4)»TC4) 

DATA S,T/ 3*0. 0,0. 5, 0.3, -0,3 ,2*0.0/ 
CALL PLOT CO. 0,0. 0,3) 

DC 1 1=1,4 

1 CALL PLOT ( SC I >,T{I),2> 

CALL PLCTC0.0,XCl)*4.0,3> 

DO 2 1=2, NX 

2 CALL PLOTC FLOAT C I ) /256.G , X( I )*4.0»21 
CALL PLOT (5 .0,0.0, 3) 

DC 3 1=1,4 

3 CALL PLOT <5 .0- S CI),TC I) ,2) 

RETURN 

ENO 
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SUBROUTINE BBPLOTIPPN) 

C — BBPLOT ASSUMES THE DATA «X" HAS THE RANGE 0.0 TO 1.0. 

COMMON /A/XI12801 ,NX»$RMS 
DIMENSION FI31) f G( 31), LABEL! 201 

f AND G CONTAIN INFORMATION USED TO PRODUCE SCALES- — 

-LABEL WILL BE USED TO LABEL THE PLOT 

DATA F, G/0.,.2* 2*0. , .1,2*0. ,.1,2*0.,. 1,2*0. ,.1,2*0. , 

1 . 1 , 2* 0 . » . 1 , 2*0 . , . 1 * 2*0. , . 1 , 2*0. , . 1 , 5*0. » 3* • 1 , 3* « 2 , 

2 3*. 3, 3*. 4, 3*. 5, 3*. 6, 3*. 7, 3*. 8, 3*. 9, 1.0/ 

CALL PLCTI0.0,8.0*XI1I,3) 

DO 1 1=2, NX 

1 CALL PL CT l F L G AT 1 1 ) / P PN » 8 • 0 *X ( 1 1 » 2 ) 

CALL PLOT (0.0, 0.0, 3) 

I FP N= P P N 
K=INX-1 ) / IPPN+l 
DO 2 J= 1, K 
8=J-1 

DC 2 L= 1, 31 

2 CALL PLOT I8+GIL ) ,-FIL ) , 2) 

XL=B+1. 0 

CALL PLOT (XL, 8.0, 2) 

DO 3 J = 1 , K 
B=J- 1 
G=XL-B 

DC 3 L = 1 ,31 

3 CALL PLOT (B-G1LI ,8.0+FIL i,2) 

CALL PLOTtO. 0,0. 0,21 
R EAD( 5 , 100 ) N 
READ! 5 , 101) LABEL 

100 FORMAT ( 15) 

101 FORMAT! 20 A4) 

SL=.14*FL0ATCN) 

T=!XL-SL)/2.0 

CALL SYMBQLIT,-.5, . 14, LABEL ,0.0, N) 

T= I XL-2. 941/2.0 

CALL' NUMBER IT, -1.0, .14, PPN,0.0»-1 ) 

CALL SYMB0LI999. 0,-1. 0,. 14,* POINTS PER INCH* «O.O t 161 
CALL PLOT (XL+1. 0,0. 0,-3) 

C MOVE PEN AND ORIGIN 1 INCH BEYOND LEFT LIMIT OF PLOT. 

RETURN 

END 
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SUBROUTINE CCCNVC 

COMMON I X , R M S D » C N A X , ICMAX,NGUTSD/A/RES( 1280 ), NS, SR MS 

COMMON /B/SP£C( 1280), R ESP C 150), NF 

RE AC (5, IOC) L COP C,N$MO, ICC PI ,N$«l , ALPHA 1, L0GP2, 

2 ALPKA2, L0CP3, ALPHA3 

100 FCRNATI4I 5,F5.2 ,2 U5,F5.2 ) } 

WRITE16,217) 

217 FORMAT ( 1H1 ,20X, 18HCECCNVCLUTICN OATA) 

WRITE {6 , 199) L OOP C » N SMO » L OOP I » NS. Ml , ALPHA 1 , L00P2 , 

1 ALPFA2* L0QP3, ALPHA3 

199 FORMAT ( 5 X6H LOOPO= I 5 /6X5HN SM 0= I 5 /5 X6HL0Q P 1 = I 5/ 

1 6X5HNSM1* 15 /4X7H ALPHA 1 = F 5.2/ 5X6H LOOP 2«*I 5/ 

2 4X7HALPHA2=F5.2/5X6HLCGP3=I5/4X7HALPHA3=F5.2/// ) 
> = C . 

CC 1 I = 1 , N F 
I F (RESP (I ) *LT*X) GC TC 1 
X= RES P { I ) 

I >=I 

1 CONTINUE 

3 IF (LCCPC. EQ.O) GO TO 5 
CG 4 I = 1 , LOOP 0 

CALL SMTHER(NSMO) 

4 WRI TE ( 6 ,200 ) I , SR PS 

200 FORMAT < / 5X 7H I LOOPO= I 5 , 5 X , 22 HRM S SMOOTH CORRECTION^, 

1 E12.5I 

5 CC 50 1=1, NS 

50 SPEC 1 1 ) =RES (,I ) 

IFtLCCPl. EQ.O) GO TO 7 
CC 6 1 = 1, LOOP 1 

CALL DCCNER I ALPHA1 ,7) 

WRITE ( 6,201 ) I » RMSD » CM AX » I OMAX » NQUTSD 

201 FCPM4T(/4X7HlLnOPl = I5,5X14HRwS DEV I AT I GN=E12. 5, 

1 5X1 6HMA X I MUM D 0 V I A T I GN = E 12 . 5 , 911 AT PCINTI5, 

2 5X2 7H NUMBER CF POINTS >1 OR <0 =15) 

CALL SMTh ERIN SMI) 

6 WRITE 16 ,202 ) I , SRNS 

202 FCPMAT</5X7HIL0CP1=I5,5X,22HR«S SMOOTH CORRECTION*, 

1 E12.5) 

7 IF{ LCCP2.S0.0) GO TC 9 
CG 8 I- 1 » LOOP 2 

CALL DCCNER ( ALPH A2 , 4 ) 

8 WPITE<6,203) I , RMSC » OMAX , I OMAX , NOUT SC 

203 FORMAT ( /4X7HI L0QP2= 15, 5X14HPMS DEV I AT IQM = E 1 2. 5, 

1 SXlEhMAXlMijv 0 5 V I A T I 0 N = E 1 2 . 5 , 9H AT PCINTI5, 

2 5X27HNUM3ER OF POINTS >1 OR <0 =15) 

9 IF(LCCP3»EQ«0) GC TO 11 

CC 10 1 = 1, LOOP 3 

CALL DCCNER ( ALPHA3 , 1 ) 

10 WPITE(6»204) I , R M S 0 i CM A X , I 0 W A X , NQUTS 0 

204 F C F M AT t /4 X 7 H J l 0 0 P 3 = 1 5 » 5 X 1 4 HR M $ 0 1 V I A T 1 0 N = E 1 2 . 5 , 

1 5X13HMAXI.MUM DEVIATION = £12*5, 9 H AT PCINTI5, 
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2 5X27HNUMBER 
11 RETURN 
END 


OF POINTS >1 OR <0 =15) 
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SUBROUTINE DCONERI P ARAM, KDELT) 

COMMON I X ,RMSD, DMAX ♦ I DMAX, NGUTSO/A/RES ( 1280 ) » NS » SRMS 

COMMON / 8/ SPEC I 1280 ) ,-RESP U50 ) ,NF 

RMSD=0. 

DMAX=0» 

NGUT S0=0 
M=NF-KDELT 
I = IX— KDELT 
12 I =1 +KDELT 
M=M+KDELT 
PRCX=0 . 

00 3 J=1 , NF 

JJ=J+I-IX 

3 PROX^PROX+RESPI J)*RESI JJ) 

C=SP EC ( n-PROX 

IF! ABSI DJ.LE. ABSIDMAX)) GO TO 77 
CMAX=D 
I DM A X = I 

77 RMSO=RMSD+D*D 
D=D*PARAM 

RES! I J-RESI I ) +D 
IFCKOELT.EQ.U GO TO 78 
KK=KDEIT-1 
DO 6 J= 1 * KK 

Q=FLGAT IKCELT-J l/FLOAT ( KDELT J*D 
RES 1 1 + J ) =RE SI I+Jl+Q 

6 PES(I-4J=RESI I-Jl+Q 

78 IFIM+KDELT. LE.NS) GO TO 12 
Z=I I— I X ) /KDELT-*- 1 
RMSD=SQRTIRMSD/Z) 

DO 7 J= IX»I 

IF (RES ( J)*RESIJ 1 .GT.RESi Jl) NOUTSD=NOUTSD+1 

7 CONTINUE 
RETURN 
END 



APPENDIX C 


The program for the deconvolution of the H 2 lines follows 
essentially the same procedure as the program for determining the 
Response function described in Appendix B. The same subroutines are 
used with the exception of GAUSSIAN and therefore only the main program 
is listed. Within the main program the base line is read, averaged and 
subtracted from each line used. After converting to a fractional 
absorption, using the value of the 100 percent absorbing point and the 
polarizer efficiency, the data is aligned and averaged as before. The 
response function is then read from cards produced by the response 
function program and the line is deconvoluted. The data is then 
converted to an absorption coefficient curve using equation (57) given 
in Section 2.5. The width at half maximum of this curve is measured 
in units of data points by the portion of the program which locates 
the 50 percent peak absorption points for the alignment of the data. 

The total area under the curve is data point units is also measured in 
100 data point intervals and totaled. 
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C THF PROGRAM TAKFS UP TC 10 RUNS OF EACH LINE F PGM 
C MAGNETIC TAPE AN C AVERAGES THEM AND OECONVQLUTES THE 
C RESPONSE FUNCTION FROM THE AVERAGE. FOR THE FIRST 
C LINE, TWO SPECIFIED BLOCKS ARE READ AND AVERAGED 
C FOR THE EASE LINE. THE INPUT CAT A IS READ AND SMOOTHED 
C ONCE. THE BASE LINE IS SUBTRACTED AND THE DATA IS 
C CONVERTED TQ A FRACTIONAL ABSORBANCE. SUCESSIVE DATA 
C IS TREATED IN THE SAME WAY AND THE LINES ARE ALIGNED, 

C AVERAGED, AND DFCONVOLUT EC . THE FINAL PROFILE IS 
C CONVERTED TO AN ABSORPTION COEFFICIENT CURVE. THE 
C HALF WIDTH AND THF AREA OF THE CURVE IS MEASURED BY 
C THE PFCGRAM. CALCC M P PLOTS ARE MADE AT INTERMEDIATE 
C STEPS TO CHECK THE PROGRESS OF THE PROGRAM. 

C DATA CARDS AS FOLLOWS: 

C CARD 1, LABEL FOR RUN ANYWHERE ON CARD 
C CARO 2, * OF LINES TO AVERAGE IN 5, 4 OF POINTS IN 
C RESPONSE FUNCTION IN 10 

C CARD 3, 10 THE 100% ABSORBING POINT IN 15, EFF THE 

C POLARIZER EFFICIENCY IN 30 

C CARD A, NUMBER Of BLOCKS READ FOR FIRST BASE BLOCK 
C {=11 IN 5, STARTING BLOCK IN 10 

C CARD 5, SAME AS 3 FOR SECOND BASE BLOCK 

C CARD 6, 8 OF BLOCKS READ FOR LINE(=10) IN 5, STARTING 

C BLOCK IN 10 

C NEXT CARD, REPEAT 3 TC 5 FOR EACH LINE 

C NEXT CARD, 4 OF CHARACTERS IN FIRST BBPLCT LABEL IN 5 

C NEXT CARD, LABEL FOR ABOVE STARTING IN ZERO 

C NEXT CARDS, RESPONSE FUNCTION 
C NEXT CARD, DECONVOLUTION DATA AS FOLLOWS: 

C U OF PASSES IN LOOP 0 IN 5, SMOOTHING INTERVAL IN 10 

C EXPRESSED IN DATA PCINT$(CCD 4 ONLY), U OF PASSES IN 

C L COP 1 IN 15, SMOOTHING INTERVAL IN 20, DECONVOLUTION 

C PARAMETER IN 25, 8 CF PASSES IN LOOP 2 IN 30, 

C DECONVOLUTION PARAMETER IN 35, 4 OF PASSES IN LOOP 

C 3 IN 40, DECONVOLUTION PARAMETER IN 45. 

C NEXT CARD, 8 OF CHARACTERS IN SECOND BBPLCT LABEL IN 5 

C NEXT CARD, LABEL FOR ABOVE STARTING IN ZERO 

C NEXT CARO, H OF CHARACTERS IN THIRD BBPLCT LABEL IN 5 

C NEXT CARD, LABEL FOR ABOVE STARTING IN ZERO 

R E A L * 4 ID 

COM MON / A / X ( 2 560 ) ,J»$RMS/t3/Y( 2 560) , RE SP ( 536 ) , NF 
DIMENSION I RU F ( 664 ) ,0(2560 ), NAME ( 20 ) , AR E A { 20 ) 

INTEGER I'M DUO) 

CALL PLOTS ( I BUF ,664,30.0, * BOYD V ’ , •RSM0TE2# 1 ) 

CALL PLOT (0.0,0. 5,-3) 

READ (5 , 723 ) NAME 

723 F C P w A T (2QA4) 

WRITE. ( 6,7 24 )NAME 

724 FORMAT (1H1,20A4) 

900 FCRMATI 1H , 10F12.5) 

132 F C F v A T ( 2 1 5 J 
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613 FCRMATl 1HC , 15X ,3HN9 = , 1 3 , 10X ,4HNBl = , I 3 ) 

725 FORMAT ( 1HO, 25H* *-*"****- **- : **********=**) 

J=2560 

CC 300 1=1, J 
800 C(I )=0.0 

REACI5, 132) ITEST.NF 
R E 6 C (5, 133) I0,5FF 
133 FORMAT ( 2 F 15.0) 

WRITE {6,502 ) 10, EFF 

502 FORMAT ( 1HC,20X ,3FIC=,F1Q.6,5X,4HEFF=,F10.6) 
W D ITE(6,725) 

C READ EASE LINE FROM TWO SEPARATE BLOCKS 
I A ,V = 1 
714 I G C = 1 

B A$E=0. 0 

245 F E A C ( 5 , 132) N 9, N 31 
WFITECt,618)NB,NBl 
CALL PEACE R{N3,N81,0tS15) 

DC 712 1=1,256 
712 E AS E = BASE + X { I ) 

IF { IGQ.EQ.2 )GQ TC 247 
IGC = 2 
CC TC 245 

247 PASE=3A$E/(512. 0*1024.0) 

WR I T E ( 6 , 3 3 3 ) B A S E 

838 FCFVAT(1H0,2OX t 5HB A S E =, F 10 . 6 ) 

C READ INPUT DATA 

BEAC(5,132)NB,M91 
J=N 3*256 

WRITE (6 ,6 19 ) !M R , N B 1 

619 FORMAT ( 1 H 0 , 1 5 X » 3 HNR = , 1 3 , 1 OX ,4HNB 1 = , I 3 ) 

WRITE (fi» 701) 

701 FCFMAK 1H0, 17HINPUT INTENSITIES) 

CALL READER (NB,NB1, 0 , 815) 

CC 938 1=1, J 

938 XII ) = X I I ) /O. 1024E + 04 
C PLCT RAW CAT A 

CALL 8 E l C C K t X » J ) 

CALL PLOT IS, 0,0. 0,-3) 

C SMCCTfH INPUT CAT A 

CALL S«THEP, (199) 

C SUBTRACT BASE LINE ANC CONVERT TO FRACTIONAL ABSORPTION 
p I C = ( I 0-8 A S E ) * E F F 
370 CC 733 1= 1, J 
788 XU ) = l X ( I ) - 6 A S E ) / RIC 
C WRITE IN PL 7 INTENSITIES 

901 WRITE (6 ,900 ) ( X ( I), 1 = 1, J, 20) 

C PLCT SMOOTHED DATA 
CALL PBLCOlUXtJI 
CALL PLC.T (6. 0,0. 0,-3 ) 

C FINC MID PC I NT 



884 GREAT=0.0 

CC 718 1=1, J 

IF (X t I ) .LT. GREAT )G0 TC 71R 
GREAT=X(I ) 

I G RE AT = I 
718 CENTIME 

HAL F=G& EAT/2. 0 
I = I C> R E A T 

768 1=1+1 

IF(XU)-HALF) 770,770,763 

770 I H A L F 2 = I 
I = I GRE AT 

769 1=1-1 

IF(X( I J-HALFJ 771,771,769 

771 IHALF 1=1 

I F { IG0.EQ.3 )GO TC S85 

I M I C ( I AM } = { IHALF1+ IHALF2 )/2 

hRI TE { 6 ,7 72 ) I GRF AT, IHALF1 , I HALF 2 , IM I D ( I AM ) 

77 2 FORMAT { 1H0 , 7HIGR F AT = , 15, 1 OX , 7H I HAL FI = , 1 5 , 10X 
1 7H I H A L F2 = , I 7 , 1 OX , 5 H I V I D= 17) 

WRITE! 6, 725 ) 

I R E DO = I M I C ( 1 ) - I M 1 0 ( TAN) 

IFClRtDC) 70,71,72 

70 CC 969 1=1, J 
K = 1-1 REDO 

IFtK.LT.l .CR.K.GT. J )G0 TO 657 
X { I ) = X ( K ) 

GO TO 569 
697 X(I)=O.OEO 
9 69 CONTINUE 
GC TC 71 
72 OC 658 1=1, J 
I = J-L + 1 
K = 1-1 RE CO 

IF1K.LT.1.QP.K.GT. J)GC TO 699 
X U ) = X { K ) 

GC TC 693 
699 X(I ) = Q . C E 0 
658 CONTINUE 
C PLCT ALIGMED DATA 

71 CALL BRLOCK (X, J ) 

CALL PLC T I-12. 0,0. 5,-3) 

C ACC LINE PROFILES 
DC 570 1=1, J 
970 C I I > = C( I ) +X( ! ) 

IF ( IAMEQ. ITEST )G0 TC 713 
I AN = I A*'+l 

C RESTART FCR NEXT LINE 
GC TC 714 

C AVERAGE LINE PROFILES 

713 CALL FLCT(18.0,-FL0AT(ITEST-2)*0.5,-3) 
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C = I TEST 
CC 971 ! — 1 * J 

971 XI )=0 { I) /C 
WRITE(6,7S2) 

782 FCPMAT ( 1H0, 21HAVERAGFD LINE PROFILE) 

WRITE! 6, 9CQX X{ I ), 1*1, J, 20) 

C PL C T AVERAGE LINE PROFILE 
CALL 8BPLCT (512*0 
C RE AC RESPONSE FUNCTION 

REAC(5,139) (RESPII) ,1 = 1, NF) 

139 FORMAT ( BE 10.4) 

WRI TE (6 ,726) 

726 FCRMAT* 1H0, 17H&ESPGNSE FUNCTION) 

WRI TF (6,900) (RES P(I) , 1 = 1, NF, 10) 

C DECCNVCUTE LINE PROFILE 
CALL OCCNVO 

C PLCT DECGNVOLUTEO LINE PROFILE 
CALL BPPLCT ( 512.0 
WRITE (6,780 ) 

780 FORMAT ( 1H 1, 25HDECCNV0LUTED LINE PROFILE) 

WRITF (6 ,900 ) ( X( I ) , I=1,J,20) 

C CALCULATE ABSORPTION COEFFICIENT 
871 CC 110 1=1, J 
110 XII) =- ( A LOG ( 1 . 0-X ( I ) ) ) 

WRIT E (6, 5 C : 5 ) 

505 FORMAT (1H1, 22HA8SCRFT ICN COEFFICIENT) 

102 KRITE(6,9C0) (X(I),I=1,J) 

C CALCULATE AREA 

DC 2C1 "=1,20 
2 01 A RE A (M) =0.0 
N = 20 
M = 1 
K=1 
L = 1 28 

WRITE<6,372) 

372 FORMAT (1H1 , 39HCALCULATEC AREAS IN HALF BLOCK SEGMENTS) 

200 DO 3C0 I = K » L 

300 AREA (M) =A PEA{M)+X( I ) 

WRI TE (fc,3C2 ) M, AREA ( N ) 

302 FCPMAT (1H , 5HAR EA ( , I 2 , 2H) * , F 15. 5 > 

N = M+ 1 
K = L 

L =K + 12S 

IF(M.EC.N+1)GC TC 202 
GO TC 200 
202 8IGA=0. 0 

CC- 400 M= l , 20 

400 BIGA=BIGA+AREA{ M ) 

WRITE ( 6, 4 C 1 ) RIGA 

401 FORMAT (IHO, 5X, 1 1HT0TAL AREA=,F15.5) 

C FIND HALF WICTF 
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I GO -3 
GC TO 8 £4 

885 10IFF=I HALF2-IHALF 1 

WRI TEC 6 ,932) I GREAT, I HAL FI, IHALF2 

93 2 FORMAT ( 1H1, 7HIGREAT*, 15 , 10X , 7HIHALF1= , I 5 , 10X, 
1 THIHALF2-, 15) 

WPITE<6,773)ID!FF 

773 FORMAT UH0,26HHALF WIDTH IN DATA POINTS=,I5J 
IFIGREAT.LT.l.OJGO TO 764 
DO 719 1*1, J 

719 XII l=XU) /GREAT 

764 WRITE 16,720) 

720 FORMAT ( 1H0, 20HRE5CALE0 ABS . COEFF.) 
WRITE(6,900)< .xm,I*L,J,20) 

C PLOT ABSORPTION COEFFICIENT 
CALL BBPLOT (512.0) 

CALL PLOTCO. 0,0. 0,999) 

15 CALL EXIT 
END 
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